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FOREWORD 


This  report  represents  one  phase  of  an  effort  aimed  at  extending  the 
techniques  of  airframe— human  pilot  systems  analysis  used  for  the  deriva¬ 
tion  of  fundamental  vehicle  handling  qualities.  Previous  investigations 
have  studied  lateral  and  longitudinal  piloting  situations  applying  single¬ 
loop  feedback  analysis  methods  and  employing  a  mathematical  model  of  the 
human  controller  as  a  servo  element.  The  multiloop  systems  analysis  tech¬ 
nique  developed  herein,  although  intended  for  manual  control  applications, 
is  equally  suitable  for  automatic  control  system  investigations. 

The  research  reported  was  sponsored  by  the  Flight  Control  Laboratory 
of  the  Aeronautical  Systems  Division  (now  part  of  the  AF  Flight  Dynamics 
Laboratory,  Research  and  Technology  Division)  under  Project  No.  8219, 

Task  No.  821905.  It  was  conducted  at  Systems  Technology,  Inc.,  under 
Contract  No.  AF  33(61  6) -8024,  with  Mr.  I.  L.  Ashkenas  and  Mr.  D.  T.  McRuer 
serving  as  principal  investigators.  The  Air  Force  project  engineer  was 
Mr.  R.  J.  Wasicko. 

The  authors  gratefully  acknowledge  technical  and  editorial  contribu¬ 
tions  made  by  Messrs.  R.  L.  Stapleford  and  R.  J.  Wasicko  and  the  careful 
work  of  the  STI  production  staff. 


ABSTRACT 


TThe  multiloop  vehicular  control  system  analysis  technique  developed 
is  designed  to  maximize  the  ransfer  to  the  multiloop  problem  of  knowl¬ 
edge  and  insights  obtained  irom  elementary  single-loop  vehicular  control 
system  analyses.  The  matrices  representing  the  closed- loop  multiloop 
system  are  expanded  in  a  special  fashion  to  forms  in  which  the  elementary 
single- loop  systems  explicitly  appear.  In  the  course  of  the  development, 
the  concept  of  the  vehicle  coupling  numerator  is  introduced  as  an  addi¬ 
tional  property  of  vehicle  response  behavior.  The  over-all  closed-loop 
system  analysis  consists  of  successive  closures  of  single,  elementary 
system  loops  and  loops  involving  vehicle  transfer  function  numerators  and 
coupling  numerators. 


The  method  is  developed  initially  using  a  general  multiloop  system. 
Two  simple,  practical  multiloop  systems  in  aircraft  control  are  then  used 
as  examples  of  the  analysis  technique;  one  illustrates  a  multiloop  situa¬ 
tion  utilizing  a  single  vehicle  control  input  and  the  other  involves 
multiloop  control  using  two  vehicle  control  inputs. 


This  technical  documentary  report  has  been  reviewed  and  is  approved. 


Chief ,  Control  Criteria  Branch 
Flight  Control  Division 
AF  Flight  Dynamics  Laboratory 
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Notational  Rules  for  Closed-Loop  Quantities 


' I  •  The  number  of  primes  present  indicates  the  number  of  loops  closed  previously 
which  affect  the  quantity  considered. 

2.  The  notation  for  the  closed- loop  factor  is  the  same  as  that  for  the  open- 

loop  factor  (plus  a  prime)  when  the  closed- loop  and  open-loop  transfer 
function  factors  have  the  same  form.  In  this  case  the  origin  of  the 
closed-loop  factor  is  always  at  hand  (e.g.,  —  qj,  ov  ov,,  q„  uu, 

Tr  Tr,  Tqj^  *  T'qj-j  >  T^  p  T^  ,  etc . )  . 

3.  When  the  closed-loop  factors  differ  in  form  from  their  open- loop  origins 
several  possibilities  exist: 

a.  For  closed-loop  factors  which  have  the  same  form  as,  and 
are  approaching,  open- loop  zeros,  the  closed-loop  factor 
notation  is  that  of  the  open-loop  zeros  (plus  a  prime) . 

For  example,  open- loop  quantities  (s  +  1 /Ts)  and  (s  +  1 /Ta2) > 
which  couple  to  form  a  quadratic  approaching  the  open-loop 
zeros  of  (s2  +  2£ro>rs  +  o^) ,  would  give  rise  to  a  closed- 
loop  factor  ordinarily  denoted  as  (s^  +  2^-^o^.s  +  . 

b .  For  closed-loop  factors  which  differ  in  form  from  both  the 
open-loop  pole  factors  from  which  they  depart  and  the  open- 
loop  zero  factors  which  they  ultimately  approach,  a  special 
notation  is  coined  which  ordinarily  reflects  the  origin  of 
the  factor.  For  example,  closed- loop  poles  which  start 
from  s  =  0  and  s  =  -1 /Tr,  then  couple  to  form  a  quadratic 
factor,  and  subsequently  decouple  to  end  finally  at  two 
real  zeros,  would  be  denoted  as  s^  +  2£rcdrs  +  odr^  in  the 
quadratic  region. 
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o.  Closed-loop  factors  which  have  no  readily  identified  origin 
or  end  point,  such  as  one  starting  at  s  =  0  and  approaching 
s  =  °°  as  gain  increases,  are  given  a  specially  coined  nota¬ 
tion,  e.g.,  l/T^. 


4.  When  the  application  of  these  rules  "by  rote  would  result  in  confusion  in 
the  local  context,  a  new  form  is  substituted  for  the  closed- loop  factor 
involved.  Primes,  however,  are  always  retained. 
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SECTION  I 


INTRODUCTION 

The  intelligent  treatment  of  linear  multiloop  systems  requires  an  analysis 
technique  which  is  comparable,  in  providing  insight  for  synthesis  activities, 
to  the  graphic  and  physically  meaningful  techniques  available  for  single-loop 
analysis.  Such  provision  of  synthesis -oriented  insights  and  physical  appreci¬ 
ation  is  central  to  the  intelligent,  interpretative  analysis  necessary  to 
achieve  near-optimum  synthesis  of  complex  systems.  Perhaps  the  most  important 
class  of  such  complex  multiloop  systems  involves  control  of  vehicles  having 
many  degrees  of  freedom.  In  these  systems  the  vehicle  is  usually  the  most 
complex  and  troublesome  dynamic  element  in  the  system;  so  much  so  as  to  ordi¬ 
narily  be  the  object  of  separate  study  and  technology.  Many  of  the  loops 
required  in  the  control  system  are  primarily  intended  to  correct  the  vehicle's 
uynamic  deficiencies.  Consequently,  the  adequacy  of  any  multiloop  analysis 
method  evolved  will  strongly  depend  on  how  easily  and  creatively  the  available 
knowledge  of  vehicle  dynamics  can  be  utilized  in  the  evolution  of  the  multiloop 
control.  Most  of  the  existing  multiloop  analysis  techniques  (see  Ref  1  for  a 
partial  summary)  essentially  ignore  this  critical  feature,  and  thereby  fail  to 
take  full  advantage  of  vehicle  stability  and  control  technology  and  the  insights 
it  can  bring  to  vehicular  control  problems. 

The  most  common  technique  for  the  pencil-and-paper  analysis  of  multiloop 
vehicular  control  systems  is  the  method  of  equivalent  stability  derivatives 
(Refs  2—6) .  Equivalent  stability  derivatives  are  especially  useful  when  com¬ 
bined  with  literal  expressions  for  the  approximate  factors  of  the  vehicle 
transfer  functions  (Refs  2  and  5) .  However,  the  method  is  restricted  practi¬ 
cally  to  ideal  (no  lag)  or  nearly  ideal  controllers,  and  it  is  most  easily 
applied  (although  not  limited)  to  control  feedbacks  which  augment  existing, 
rather  than  create  new,  derivatives.  Within  these  limits  it  is  an  extremely 
useful  analysis  and  synthesis  tool. 


Original  manuscript  submitted  December  1962;  revised  manuscript  released  by 
authors  January  1964  for  publication  as  an  ASD  Technical  Documentary  Report. 
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By  far  the  most  common  means  for  multiloop  system  study  and  synthesis  i-s 
repetitive  analysis  using  the  analog  computer.  With  the  aid  of  this  tool  and 
a  background  comprising  knowledge  of  the  vehicle  equations  of  motion,  an  appre¬ 
ciation  of  the  gross  single-loop  effects  of  probable  feedbacks,  approximate 
factors  for  the  vehicle,  past  computer  solutions,  etc.,  the  analyst  can  rela¬ 
tively  rapidly  arrive  at  a  set  of  loop  closures  which  define  a  "good"  system. 

A  "good"  system  in  this  respect  will  exhibit  fast,  well-damped,  accurate 
responses  to  all  representative  commands  and  will  act  similarly  in  suppressing 
disturbances  (i.e.,  responses  to  commands  and  disturbances  will  be  akin  to 
those  of  a  well-damped,  low-order  system) .  The  "good"  system  will  also  be 
relatively  insensitive  to  changes  from  nominal  values  in  the  vehicle  or  con¬ 
troller  characteristics,  will  tend  to  be  unaffected  by  the  introduction  of 
gmall  parasitic  nonlinearities,  etc.  With  the  computer  there  is  no  essential 
difficulty  in  treating  controller  lags  or,  for  that  matter,  nonlinearities. 

Exclusive  reliance  on  the  analog  computer  does  have  some  deficiencies.  For 
instance,  some  insight  into  the  over-all  system  is  irretrievably  lost  because 
of  the  dominance  of  only  particular  modes  in  the  time-histories.  Also,  modes 
which  may  be  of  great  importance  when  the  conditions  are  changed  slightly  are 
stcpp ressed,  sensitivities  are  difficult  to  evaluate,  gross  trends  and  grand 
simplifications  are  harder  to  come  by  than  with  some  analytical  schemes,  etc. 
Finally,  elements  described  only  in  frequency  response  terms,  such  as  experi¬ 
mentally  measured  subsystem  describing  functions  and  human  pilot  descriptions, 
cannot  be  used  directly  in  computer  operations.  Thus  insight  is  constricted 
and  initiative  stifled,  as  always  happens  when  only  a  single  approach  to  a 
problem  is  utilized. 

The  properties  desired  of  a  multiloop  analysis  technique,  in  the  light  of 
the  above  renarks,  should  include: 

1 .  Analytical  formulations  which  show,  as  separate  entities, 
vehicle-alone  and  controller-alone  characteristics 
expressed  in  conventional  and  well-understood  terms — 
thereby  providing  a  close  tie  to  the  individual  elements 
and  the  physical  problem. 

2.  Analytical  operations  which  can  be  performed  using  the 
more  efficient  of  the  classical  graphical  techniques  of 
servoanalysis — thereby  enhancing  transfer  of  skill  and 
intuition. 
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3°  Analysis  sequences  and  procedures  which  are  highly  respon¬ 
sive  to  physical  insights  and  intuition;  and  which  lead, 
when  used  by  a  skilled  practitioner,  to  "good"  systems  with 
a  minimum  of  iteration. 

4.  Problem  solution  presentations,  and  results,  which  are 

supplementary  as  well  as  equivalent  to  the  results  obtained 
using  the  analog  computer. 

This  report  develops  and  explains  a  multiloop  analysis  technique  which  is  intended 
to  satisfy  most  of  the  wants  discussed  above. 

Like  almost  all  "new"  methods  the  technique  advanced  here  has  some  precedents 
which  were  significant  in  its  evolution.  Various  steps  in  the  technique  to  be 
explained  were  evolved  over  a  fairly  long  period  of  time,  starting  about  1957  at 
Systems  Technology,  Inc.  An  analysis  technique  used  in  Ref  7  played  an  important 
role  in  suggesting  some  aspects  of  the  development.  Also,  an  analysis  procedure 
having  some  features  similar  to  that  presented  here  was  used  by  Mataga  in  Ref  8. 

In  many  respects  the  unified  servoanalysis  procedure,  as  reported  in  Ref  9, 
and  the  sensitivity  and  modal  response  analysis  techniques  reported  in  Ref  10, 
are  the  analytical  companions  to  the  multiloop  analysis  techniques  reported  here. 
This  trio  of  reports  ir  intended  to  cover,  with  a  unified  eclectic  point  of  view 
regarding  methods,  most  of  the  significant  analysis  problems  of  linear  servo 
theory. 

The  body  of  the  report  is  presented  in  three  sections  .  The  major  analytical 
effort  appeal's  in  Section  II  for  a  sequence  of  generalized  systems  .  Generalized 
equations  of  motion  are  used  as  a  starting  point,  then  formulated  as  matrix  equa¬ 
tions  to  simplify  much  of  the  analytical  development.  Certain  key  observations 
on  the  desired  types  of  matrix  expansions  are  then  made,  after  which  the  analyti¬ 
cal  development  is  largely  straightforward  algebra.  At  the  end  of  this  section 
some  practical  problems  of  application  deriving  from  physical  considerations  are 
listed  and  discussed  in  general.  The  last  two  sections  of  the  report  are  intended 
to  elucidate  by  example  some  of  the  aspects  of  the  method  which  derive  predomi¬ 
nantly  from  these  physical  considerations.  These  sections  start  with  the 
generalized  treatment,  and  specialize  it  for  two  concrete  examples.  The 
discussion  of  Section  III  treats  a  longitudinal  example  for  an  altitude  control 
system.  This  particular  system  is  the  simplest  possible  multiloop  system 
insofar  as  its  basic  loop  structure  is  concerned.  The  lateral  example,  covered 
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in  Section  IV,  is  considerably  more  complex  and  serves  as  a  simplified  prototype 
for  most  multiloop  vehicular  control  problems. 


4 


SECTION  II 


ANALYSIS  OF  GENERALIZED  MULTILOOP  SYSTEMS 

In  this  section  the  essential  features  of  a  multiloop  analysis  technique 
will  be  developed.  A  generalized  notation  for  vehicle  and  controller  variables 
is  introduced  early  and  used  throughout  the  developments.  Matrix  formulations 
are  appropriate  for  multiloop  problems,  and  could  be  used  from  the  outset;  but, 
to  make  the  developments  easier  to  follow  for  the  reader  unfamiliar  with 
matrices,  exemplary  equations  of  motion  are  intermixed  with  matrix  generaliza¬ 
tions.  In  most  of  the  developments  an  inductive  approach  is  used  wherein 
systems  of  limited  complexity  are  used  to  formulate  equations  which  are  then 
expressed  in  matrix  form.  The  matrix  equations  so  obtained  are  both  a  short¬ 
hand  for  the  limited  conplexity  system  equations  and,  viewed  more  broadly,  the 
appropriate  equations  for  far  more  complex  systems.  Thus  an  attempt  is  made  to 
satisfy  the  often  conflicting  desires  for  concreteness  and  generality. 

Several  types  of  systems  appear  in  the  course  of  the  development.  These 
differ  primarily  in  their  level  of  complexity,  which  has  been  selected  to  be 
just  sufficient  to  illustrate  the  local  points  being  made.  The  first  system 
considered  (in  articles  A  —  C)  is  the  most  complex;  it  is  used  to  illustrate 
the  generalized  notation,  matrix  formulation,  and  closed-loop  system  character¬ 
istic  equation  development.  When  closed-loop  transfer  function  numerator  terms 
become  the  subject  of  detailed  examination,  a  somewhat  simpler  system  is  intro¬ 
duced  for  purposes  of  clarity.  Finally,  in  article  D,  a  still  simpler  system 
is  used  to  explain  the  steps  involved  in  the  final  phases  of  the  analysis 
process.  The  two  simpler  systems  are  special  cases  of  the  more  complex  one, 
which  is  introduced  below. 

A.  GENERALIZED  SYSTEM 

The  generalized  vehicular  control  system  to  be  analyzed  is  shown  in  the 
block  diagram  of  Fig.  1 .  Despite  its  foreboding  appearance,  the  system  is 
relatively  simple,  being  just  complex  enough  to  represent  most  flight  control 
situations  and  to  allow  easy  inductive  generalizations.  The  system  consists  of 
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a  vehicle  plus  control  equipment  comprising  sensing,  equalizing,  and  actuating 
elements.  The  vehicle  has  three  independent  degrees  of  freedom,  and  is  subject 
to  control  forces  and  moments  applied  by  two  control  deflections  and  two  exter¬ 
nal  disturbances.  The  control  deflections  are  functions  of  command  inputs, 
feedbacks  from  the  three  degrees  of  freedom,  and  a  possible  fourth  feedback 
from  an  auxiliary  variable  which  is  a  function  of  the  independent  degrees  of 
freedom.  Additional  auxiliary  feedbacks  may  also  be  present,  but  the  single 
one  shown  will  illustrate  the  analytical  process  relating  to  such  quantities. 
Generalization  by  induction  tc  include  added  auxiliary  feedbacks  will  later  be 
seen  to  be  straightforward. 

The  Fig.  1  block  diagram  can  be  simplified  to  one  having  only  unity  feedback 
loops  by  reinterpreting  the  command,  Jj ,  and  forward  loop,  G-j_  j ,  transfer  blocks. 
This  could  be  done  in  the  q2  loop,  for  example,  by  replacing  J2(s)  with 
J2(s)  =  J2(s)/h2(s),  and  replacing  G^ 2(s)  and  G22(s)  by  G12(s)  =  Gl2(s)H2(s)  and 
G22(s)  =  G22(s)H2(s) ,  respectively.  The  unity  feedback  loops  block  diagram  of 
Fig.  2  results  when  similar  steps  are  taken  throughout.  This  block  diagram  is 
still  somewhat  overcomplicated  in  that  command  inputs  are  shown  for  all  feedback 
loops.  Actual  commands  ordinarily  exist  for  only  one  loop,  or  possibly  two. 

The  analytical  operations  involved  in  the  above  reduction  to  unity  feedbacks, 
and  many  other  operations  to  follow,  are  based  on  the  assumption  that  system 
elements  can  be  described  by  linear,  constant-coefficient,  differential  equations. 
Accordingly,  all  motion  and  transfer  quantities  in  Figs.  1  and  2  are  shown  as 
functions  of  s,  the  Laplace  transform  complex  variable.  The  assumption  of 
linearity  also  facilitates  the  manipulation  of  system  equations  and  allows  the 
principle  of  superposition  to  be  used.  Superposition  provides  a  great  conveni¬ 
ence  because  the  generalized  closed- loop  transfer  functions  need  be  developed 
for  only  one  command  input.  Closed-loop  transfer  functions  for  other  commands 
can  then  be  found  by  proper  juggling  of  subscripts  and  sui table  interchange  of 
elements  within  particular  matrices. 

To  go  along  with  the  economy  of  notation  provided  by  the  use  of  macrix 
methods,  the  functional  dependence  of  most  quantities  on  s  is  not  indicated 
hereafter.  For  example,  the  vehicle  equations  of  motion  for  the  three  inde¬ 
pendent  degrees  of  freedom  are  written  as 
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(1) 


a1  1  q^  +  ai  292 

+  a,  3q3  =  F^ 

+  F12B2  +  E11r)1 

+ 

E1  2r)2 

a2i  9-j  +  a2p92 

+  a2393  =  F2181 

+  F2252  +  Eg1  r)1 

+ 

EggTJg 

a 31  9-|  +  a 3292 

+  a539.3  =  f315-| 

+  F^qSq  +  E?1  r)1 

+ 

E32rl2 

where  the  a^j's  are  functions  of  s  and  vehicle  characteristics  such  as  stability 
derivatives  (Ref  11).  Then,  the  vehicle  transfer  functions  for  the  responses  of 
the  independent  outputs  to  control  deflections  are 


&1 

fl_ 

s2 

fa 

5i 


V<H 


b2 


A 


fi 

5J 


where  A  is  the  characteristic  determinant  of  the  vehicle  equations. 


(2) 


a1  1 

a1  2 

a1  3 

a21 

a22 

a23 

a3l 

a32 

a33 

(3) 


and  the  numerator,  Nqis_,  is  obtained  by  replacing  the  column  of  qi  coefficients 
in  A  by  the  column  of  5j  coefficients  given  on  the  right  side  of  Eq  1  . 

Additional  shorthand  defining  the  system  and  the  closed-loop  transfer  func¬ 
tions  is  also  needed.  The  information  required  to  specify  a  particular  multiloop 
system  is 
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The  command  or  input  to  the  system 

The  vehicle  output  to  be  controlled 

The  control  deflections  to  be  used 

The  error  signals  (or  feedbacks)  that 
activate  the  controllers 

This  information,  using  a  particular  system  as  an  example,  will  be  designated 
as  follows : 

q^82  ;  V  ('0 


Here  the  feedback  signals  to  52  are  and  q,  ,  and  to  6-j  are  q^  and  q^.  The 
system  input  is  a  command,  q-|  ,  which  also  implies  the  directly  controlled 
output.  Note  that  the  q^ — feedback  is  implied  by  the  q-j  ,  so  a  separate 
call-out  is  redundant.  The  closed-loop  transfer  function  for  this  system  is 
designated  by 


V  1 &1 

V  q4—"~62 


and  the  appropriate  controller  equations  are 

81  =  g12(-cJ2)  + 

(  8) 

82  =  G21  )  +  G24(-cV  +gig21<1ic 

A  disturbance,  such  as  ,  can  also  appear  as  a  system  input. 

Other  simplifying  notation  pertinent  to  specific  steps  involved  in  the 
development  of  closed-loop  system  transfer  functions  will  be  presented  later 
as  required. 

B.  FORMULATION  OF  THE  SYSTEM  EQUATIONS 

The  complete  equations  of  motion  for  the  vehicle  including  the  auxiliary 
output  motion  quantity,  q^,  can  be  considered  in  two  ways.  The  first  uses  the 
three  degrees  of  freedom  set,  Eq  1  ,  as  the  vehicle  equations  and  then  obtains 
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auxiliary  outputs  as  linear  combinations  of  the  three  independent  outputs.  For 
example,  for  the  linear  combination  would  be  q^  =  a4i  ^  +  a42^2  +  ^3^3* 

After  substituting  this  relationship  into  the  control  deflection  equations,  the 
terms  -G-|  ^a^  q-|  ,  -Gijj.ajj.gqg,  and  -G-|  ^a^^q^  will  appear  in  the  equation  for  5-| 
and  the  terms  -Ggjj.aij.iqi,  -Ggjj.aij.2q2>  and  -Ggjj.ajj.3q3  will  appear  in  the  equation 
for  5g.  When  combined  with  the  terms  representing  controller  action  on  the 
independent  outputs,  equivalent  controller  transfer  functions  are  formed.  For 
example,  GiJjajjq  +  Gi  j  ,  GiJjajjg  +  Gig,  and  Gi  Jj.ajj.j  +  Gi  3  replace  Gi  1  ,  Gig,  and  Gi  3 
as  the  forward  loop  transfer  functions  relating  q^ ,  qg^,  and  q^-  ,  respectively, 
to  5i .  The  q^  command  function  can  be  assigned  similarly  to  the  commands  for 
the  independent  vehicle  outputs.  Thus,  when  this  approach  is  used  the  auxiliary 
output  motion  quantity,  q^,  effectively  disappears  into  the  controller  portion 
of  the  system. 

The  second  method  of  including  an  auxiliary  output  motion  quantity  is  to 
treat  it  in  the  same  way  as  the  independent  vehicle  outputs.  With  this  scheme 
the  auxiliary  output  is  viewed  in  association  with  the  vehicle  portion  of  the 
system.  Confining  attention  to  the  exemplary  case,  the  three  degrees  of  freedom 
set,  Eq  1  ,  is  combined  with  the  equation  for  q^.  The  latter  is,  of  course, 
redundant  as  a  system  equation.  Thus, 


a11q1 

+  al2qg 

+ 

=  FiiBl 

+  F12B2 

+  E1 1  ’’•l 

+  F-|  2^2 

a2lq1 

a22^2 

+ 

a23q3 

=  F21B1 

+  F22B2 

+  ^1 

+  "^22^2 

(7) 

a31  al 

+  a32q2 

+ 

a33q3 

=  F3i5i 

+  F5262 

+  E31TI1 

+  E32T12 

'a4iql 

-  ajjp'Jfi 

- 

a43q3  +  q4 

=  0 

or  in  matrix  notation  (Ref  1 2) , 

[a]  DO  =  [F][5]  +  [s]  W  (8) 


ai  1 

ai2 

ai3 

0 

a2l 

a22 

a23 

0 

a3l 

** 

a33 

0 

"a4i 

"a42 

'a43 

1 

1 1 


[aj 


(9) 


(12) 


0  3) 


OM 


With  four  simultaneous  commands  to  each  controller,  the  transformed 
equations  of  motion  of  the  controllers  for  the  system  of  Fig*  2  are 


51  -  G1 1  (Ji*i_  -  )  +  Gi  2(J2q2c  "  q2)  +  G1  3^J3q‘3c  “  q3^  +  G11^JV34C  " 

(15) 

52  =  G2i  (J^c  -  Q-i  )  +  g22  ( J 2q2c  "  q2^  +  G23(J3<13C  "  q3^  +  G24^  JV14C  " 


12 


or  in  matrix  notation, 

DO  =  [<0W 

=  [G]  H  [Qj  -  [G]  W 

=  [Gel  &c]  -  [G]  [q]  (16) 


where 


s 

II 

- 1 

_p 

G12 

Gl  3 

°"1 

(17) 

LG21 

G22 

G23 

°2lJ 

W  =  [Jl 

J2 

3 

P 

1 _ 1 

(18) 

G11  J1 

G12J2 

G1  3J3 

Gl4J4 

(19) 

= 

_G21J1 

^22*^2 

g23j3 

Ggl^Jk 

With  the  formulation  of  the  matrix  equations,  the  system  block  diagram  of 
Fig.  2  can  be  replaced  by  the  matrix  block  diagram  shown  in  Fig.  3-  This  is  as 
deceptive  in  its  simplicity  as  Fig.  1  is  for  apparent  complexity.  The  system 
equations  of  motion  can  be  found  from  either  Fig.  3  or  Eq  8  and  16.  Thus, 
substituting  Eq  16  into  Eq  8,  transposing  and  collecting  like  terms,  leads  to 

[GO  +  [p]  [g]]  00  =  0]  [Gj  [qj  +  [e]  DJ  (21 ) 

and  after  premultiplying  by  the  inverse  of  [W  +  M  tel]  ,  the  explicit  expression 
for  M  becomes 


13 


lb- 


Figure  3-  Matrix  Elock  Diagram  for  General  Control  System 


(22) 


This  is  the  formal  solution  for  the  closed-loop  system.  However,  much  remains 
to  be  done  before  it  has  any  concrete  value  in  analysis  or  synthesis  procedures 
using  servo  methods. 


C.  REDUCTION  OF  THE  MATRIX  EQUATIONS 

From  Eq  9,  11,  and  17,  after  performing  the  indicated  matrix  multiplication 
and  addition. 


[a]  +  [F][G] 


*1  1 

+  F-,  i  G-|  i 

+  Fl2G2i 


*21 

+  F21G11 
+  f22g21 


l3i 

+  F51Gn 

+  Fj2G21 


*41 


i1 2 

+  Fi  i  G1  2 
+  F12G22 


*22 

+  f21g12 
+  F22G22 


"32 

+  f51g12 
+  f52g22 


a42 


l1  3 

+  FnG13 


+  F- 


1  2g23 


*23 


r21u1  3 
+  f22g25 


l33 

+  F51  Gt  5 
+  F52G25 


F1  i  G1 4 
+  f12g24 


F21  G1  4 

+  f22g24 


f31  gi  4 

+  Fj2G24 


-a4j 


(23) 


When  set  equal  to  zero,  the  determinant  of  Eq  25  is  the  closed-loop  system  char¬ 
acteristic  equation;  for  reasonable  controller  dynamics  (e.g.,  G±j  a  second-  or 
third-order  system),  and  a  normal  vehicle  (e.g.,  a-j_j  first-  or  second-order  in  s)  , 
it  can  be  a  fearsome  thing  indeed.  The  routine  application  of  digital  computing 
methods  can  be  used  to  factor  the  characteristic  equation  when  all  elements  are 
specified  numerically.  However,  in  a  synthesis  procedure,  the  G's  are  unknown 
within  a  wide  range  of  possible  variation.  Computing  routines  which  cover  such 
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ranges  in  increments  have  been  used,  but  the  mass  of  data  to  be  digested  can  be 
immense.  Physical  appreciation  and  the  transfer  to  the  multiloop  problem  of 
understanding  gained  with  simpler  systems  is  almost  nil.  The  net  result,  more 
often  than  not,  is  to  forget  direct  analysis  except  possibly  as  a  check  on  the 
final  results  obtained  after  using  an  analog  computer  to  narrowly  delimit  the 
ranges  of  suitable  controller  functions.  In  any  event,  pencil-and-paper  methods 
of  conventional  servoanalysis  are  not  well  suited  to  the  problem  presented  by 
Eq  21  through  23  as  they  stand. 

Four  considerations  which  relate  to  a  change  in  the  detailed  form  of  the 
problem  to  one  better  suited  for  servoanalysis  procedures  are: 

1  .  The  vehicle  is  ordinarily  the  most  complex  element  in  the 
system  and  its  dynamics  tend  to  be  dominant  - 

2.  The  study  of  vehicle-alone  dynamics  often  has  separate 
status  as  an  engineering  field  (e.g. ,  aeronautical  stability 
and  control) .  Consequently,  the  vehicle  dynamics  are  fairly 
well  understood. 

3.  Many  of  the  feedback  loops  involved  in  multiloop  systems 
exist  primarily  to  correct  vehicle  dynamic  deficiencies 
(i.e„,  parallel  equalization)  or  to  suppress  particular 
types  of  disturbances.  The  type  and  general  form  of  most 
such  feedbacks  can  be  derived  readily  from  a  knowledge  of 
the  vehicle-alone  dynamics. 

Conventional  servoanalysis  methods  are  ideally  suited  to 
close  single  loops,  e.g.,  to  find  G(s)/[l  +G(s)]  given  G(s). 

A  fruitful  approach  to  the  multiloop  problem  will  use  these  factors  to  advantage. 
The  scheme  developed  here  does  so  by  (l )  expanding  the  determinant  of  Eq  23 
in  such  a  way  that  the  vehicle  characteristic  determinant.  A,  and  the  vehicle 
transfer  function  numerators,  Wq  ,  appear  explicitly  in  the  closed-loop  system 
characteristic  equation,  and  (2)  treating  the  resulting  expressions  as  a  series 
of  "equivalent"  single-loop  servo  systems.  The  first  step  provides  vehicle- 
alone  dynamics  as  recognizable  separate  entities  in  the  system  equations  instead 
of  as  elements  in  a  hodgepodge  of  vehicle  and  control  terms.  This  allows  direct 
application  of  vehicle-alone  knowledge  and  understanding  to  multiloop  control 
problems.  The  second  step  permits  many  of  the  fruitful  insights  and  techniques 
of  single-loop  servoanalysis  to  be  applied  to  multiloop  situations. 

Expansion  of  the  determinant  of  Eq  23  in  such  a  way  as  to  retain  vehicle-alone 
transter  function  elements  as  separate  entities  is  accomplished  in  Appendix  A. 
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The  result  is: 


a1 1 
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+  F-j 
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21  G1 1 


22G21 


a22 

+  f2ig12 
+  ^22^22 


a23 

+  f21  g1  5 

+  f22g25 


F21  g1  k 

+  F22G2l, 


a31 

+  F^Gn 
+  F3gG2l 


a32 

+  f31  G,  g 
+  f32g22 


a  33  F3i  gi  4 

+  F31  G-]  3  +  FjgGgi,. 

+  f52G23 


_a4l 


-  a  *  2  2  Oji%l6j 


"a42  _a43 

+  L  £  G1iG2kN^g 

i=1  k=1  1 


(24) 


A  and  No.  are  the  characteristic  function  of  the  vehicle  and  the  numerator  of 
the  Qigi  transfer  function,  respectively.  They  are  tabulated  in  terms  of  airframe 
parameters  for  the  usual  degrees  of  freedom  and  control  deflections  involved  in 
aircraft  control  in  Refs  2,  and  11.  In  the  notation  presented  so  far  in  this 
report. 
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F2j 

%  = 

a31 

a  32 

a33 

F*. 

"a4i 

~ak2 

-ak3 

0 

(25) 


(26) 


Nq^  's  for  the  independent  outputs  (i  =  1,2,5)  will  simplify  to  5*5  determinants. 

J  ‘li'lk 

Terms  of  the  form  Ng1 gg  in  Eq  24  are  called  coupling  numerators .  They  are 
found  by  replacing  the  ith  and  kth  columns  of  Eq  9  "by  the  first  and  second 
columns,  respectively,  of  Eq  1 1  .  As  a  particular  example  using  independent 
outputs, 
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*3*2 

51  52 
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F11 

0 

a21 

F22 

F21 

0 

a  31 

f32 
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0 

_a4l 
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0 
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a11 

F1  2 

F11 

a21 

■  F22 

F21 

a51 

F52 
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(27) 


This  determinant  is  recognized  as  the  characteristic  determinant  with  the  q^ 
and  q2  column  terms  replaced  by  5^  and  52  control  effectiveness  terms,  respec¬ 
tively.  The  awkward  but  highly  descriptive  symbol  ^-n'ten‘^e^  "bo  suggest 
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this  replacement*  It  is  apparent  that  Ng  has  no  meaning  when  i  =  k*  By 
arbitrarily  defining  this  to  be  equal  to  zero,  the  "i  £  k"  note  on  Eq  24  can 
be  removed.  The  coupling  numerators  have  other  properties  of  interest  such 
as 


%l6i 

liik 

=  %262  =  0 

(28) 

NS-,52 

Nqiqk  %% 

-  -  wB251 

(29) 

Jk.%. 

\&2 

=  7|Un.  Na,  -  Wa.  Nai  \ 

A  ^  cixBl  Gkgg  li6g  q.kBl  ) 

(30) 

Common  coupling  numerators  for  aircraft  control  are  presented  in  Appendix  B. 

While  the  system  characteristic  equation,  A.^,  is  the  denominator  for  all 
closed-loop  transfer  functions,  regardless  of  command  inputs,  the  numerators 
depend  on  the  particular  command  input.  The  outputs  due  to  command  inputs  can 
be  obtained  from  Eq  22.  With  disturbances  zero, 


DO  -  [M  ♦  MCo]]’'W[oJ[qJ  (5 

where  [<J  ,  [f]  ,  and  &J  are  given  in  Eq  10,  11,  and  20,  respectively.  The 
matrix  [M  *  MM]-  ^  is,  after  inversion,  expressible  as 


[[a]  +  [FjtGj]-1 
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A3i 
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*1  2 
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a32 

*42 

*1  3 

Agj 

a33 
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*1 4 

Ag^ 

a34 

*44^ 

^sys 


(32) 


where  the  numerator  is  the  transpose  of  the  matrix  made  up  of  the  cofactors  of 
Eq  2k .  Thus , 
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*11  = 
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a33 

+  f51g15 
+  f32g23 
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F21Gl4 

+  F22G2,4. 


f31gi4 

+  F32&24 


(33) 


is  Jhe  cofactor  of  (an  +  FnGn  +  F-|  2G21  )  in  the  determinant  of  Eq  24.  Similarly, 
A]  2  is  the  cofactor  of  (a-|2  +  F^G^  +  Fi2G22)  in  the  same  determinant.  The  other 
matrix  involved  in  Eq  31  is  03  [Gj  CoJ  •  For  the  typical  case  where  [Gj  =  [G]  , 
the  matrix  multiplication  of  03  CG  3  Ocl  results  in 


[F][G3[qc]  = 
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Using  only  q-]  as  an  example  (i.e.,  q2  =  G3  =  %  =  0) ,  the  closed-loop  system 


transfer  function  is 


%  1  ^F1  1 G1  1  +  Fi2G2l^  +  ^21^F21G11  +  F22G2p  +  A31^F31G11  +  F32G21  j. 


*1, 


Asys 


(35) 


Closed-loop  transfer  functions  for  q2/q^,  etc.,  or  transfer  functions  for 
disturbance  inputs,  such  as  q.  /q-|  ,  are  obtained  in  a  similar  way. 
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Although  the  development  can  be  continued  along  the  above  lines,  the  details 
are  so  tedious  that  the  main  stream  of  the  argument  is  obscured.  Therefore,  a 
considerably  simpler  system,  shown  in  Fig.  4,  will  be  used.  (The  matrix  block 
diagram  of  Fig.  3  still  applies.)  Even  this  system  is  more  complex  than  most 
flight  control  systems,  so  still  further  simplifications  will  later  be  made. 


Figure  4.  Block  Diagram  of  Multiloop  System 
<lic  >  ~ ^ 


In  Fig.  4,  all  of  the  command  and  disturbance  inputs  are  zero  except  <T]C, 
and  there  is  only  one  feedback  to  the  controlled  variable.  The  remaining  feed¬ 
backs  are  intended  only  to  modify  the  basic  vehicle  characteristics.  The 
closed-loop  transfer  function  of  the  primary  controlled  variable,  ,  is,  from 
Eq.  35, 
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or,  written  as  determinants 
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The  denominator  can  tie  most  readily  obtained  by  specializing  Eq  24;  that  is. 


^sys  =  A  +  G1  2n02  +  G1 4^04  +  G21  Nq1  +  G24NQ4 


i81  ’  '  'H81  *'B2  ^82 

O2I4  q4qi 

+  G-|  2g21  52  +  g12g24NB752  +  Gl4G2lN5l&2 


(38) 


In  the  numerator  and  denominator  of  Eq  37,  the  cofactors  of  terms  containing 
G2i  are  identical.  Thus,  q-j  / q-j  can  be  obtained  directly  from  the  system 
characteristic  as  given  in  Eq  24  or  38.  The  closed-loop  transfer  function 
q-j  /q-j  is  then 

q4q1 ' 
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jh 
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G2l(Nq1g  +  Gl2%f52  +  Gl4N5^52') 
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(39) 


Equation  39  contains  vehicle  transfer  function  numerators  and  denominators  as 
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separate  entities.  By  dividing  the  numerator  and  denominator  of  Eq  39  by  A, 

Q.1  /‘lie  can  be  expressed  in  terms  including  vehicle  transfer  functions  themselves 
that  is. 
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where  Qo0  ,  Qk0  ,  ft.  ,  and  Qj.  are  vehicle  transfer  functions. 
c5i  ^5-|  1  b2  u>2 

The  bracketed  term  in  the  denominator  of  either  Eq  39  or  40  is  seen  to  be 
identical  to  the  numerator  terms.  This  makes  it  easy  to  recognize  that  the 
open-loop  transfer  function  q^/q-j^  is  just 
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or,  dividing  both  numerator  and  denominator  by  A, 
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Although  not  directly  apparent  because  of  their  complexity,  Eq  39  through 
42  are  of  a  form  to  permit  knowledge  of  vehicle-alone  characteristics  to  be 
used  with  maximum  benefit  in  system  synthesis  procedures  based  on  conventional 
servoanalysis  methods.  All  of  the  terms  contain  vehicle-alone  transfer  functions 
or  ratios  of  coupling  and  conventional  numerators;  and  only  controller  transfer 
functions  (the  G's)  multiply  such  terms.  Also,  the  total  transfer  function 
forms  of  Eq  39  through  42  can  be  generated  readily  by  successive  loop  closures 
(a  mere  seven  in  the  case  of  Eq  39  40!)  using  these  elemental  vehicle  and 

controller  characteristics.  Remarkably,  it  is  possible,  as  will  be  illustrated 
by  subsequent  examples,  to  rapidly  and  effectively  perform  these  operations  in 
a  fashion  leading  to  a  good  closed-loop  system  while  retaining  physical  apprecia¬ 
tion  and  developing  insights  throughout  the  process. 

Indirectly  controlled  output  transfer  functions,  qn/qic  (n  /  i),  are  also 
required  for  some  purposes .  Such  transfer  functions  as 
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^sys 


(43) 


can  be  determined  from  an  expansion  of  Eq  22.  The  characteristic  function, 

Agyg,  has  already  been  given  in  symbolic  terms  for  the  general  system  in  Eq  24 
and  for  the  simplified  system  of  Eig.  4  in  Eq  38,  so  expressions  for  Nn  alone 
are  required.  These  are  easily  obtained  from  Nj  by  replacing  q^  subscripts 
and  superscripts  by  qn's  in  the  transfer  function  numerator.  For  example, 
using  9i/q.|c  as  given  by  Eq  39,  the  closed-loop  transfer  function  q^/q^  becomes 
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This  indirectly  controlled  output  transfer  function  can  also  he  obtained  by 

,  ,  /  9k<li  \ 

multiplying  and  dividing  Eq  44  by  G21  ^  +  G-|  2Ng.)  gg  +  G-|  4%^  gg  J  and  using 


Eq  39.  The  result  of  these  operations  is 


where 


&! 

%~~b2 


%-~S2 


V  +  G1  5p 


929,  949, 

-  ~  +  Gl4%5 


1  2 


v9! 


52 


+  G-i  oN, 


929i 


949i 


2W6152  +  Gl4%1S2 


(4-5) 


(46) 


An  alternate  procedure  for  formulating  the  system  closed-loop  transfer 
function  numerator  and  denominator  expansions  is  presented  in  Appendix  C. 

D.  GENERAL  ANALYSIS  PROCEDURE 

The  developments  above  provide  general  formulations  of  the  multiloop  problem 
in  matrix  form,  and  illustrate  a  special  expansion  for  the  matrix  equations.  The 
results  of  the  special  expansion  are  various  closed-loop  transfer  functions 
expressed  in  terms  of  elemental  vehicle  and  controller  characteristics.  The 
final  step  is  the  application  of  conventional  servoanalysis  techniques  to  the 
transfer  functions  as  formulated  above.  This  phase  will  be  discussed  here  in 
general  terms,  and  the  specifics  will  be  illustrated  in  succeeding  sections  using 
particular  systems. 

To  make  the  discussion  as  simple  and  concrete  as  possible  while  retaining 
the  generalized  notation,  consider  the  system  given  by  the  block  diagram  of 
Fig.  5.  This  is  the  simplest  possible  multiloop  system  using  multiple  control 
deflections.  It  is  just  complex  enough  to  exhibit,  on  a  rudimentary  scale, 
the  major  features  of  far  more  complicated  systems.  Further  simplification, 
to  a  system  with  one  active  control  deflection  but  retaining  two  output  motion 
feedbacks,  results  in  a  multiloop  system  which  is  too  limited  to  illustrate 
all  the  features  of  interest . 
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Vehicle 

Dynamics 


Figure  5 •  Block  Diagram  of  Multiloop  System 
q.ic— *“52i  Gg— — &i 


Specialization  of  Eq  39  and  40  to  the  system  of  Fig.  5  (G^  =  G24  -  0)  results 
in  the  closed- loop  transfer  function  forms  given  below. 
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The  open-loop  transfer  function  for  the  command  loop,  with  the  other  (q2)  loop 
closed,  is  a  special  case  of  Eq  4l  and  42 ;  that  is. 
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The  vehicle  characteristics  for  this  system  are  given  by: 


a1 1 

al  2 

al  3 

A  = 

a21 

a22 

a23 

a3l 

a32 

a33 

F12 

al  2 

ai  3 

F22 

a22 

a23 

!Jq 

ql62 

F32 

a32 

a33 

A 

A 

al  1 

Fli 

al  3 

a21 

F21 

a23 

X 

a3l 

f31 

a33 

A 

A 

F1  2 

F11 

ai  5 

*2*1 

Mr  r 

51S2 

f22 

F21 

a23 

F32 

F3i 

a33 

(51 ) 


(52) 


(53) 


(5t) 


27 


An  "equivalent  block  diagram"  for  Eq  kj  through  50  is  given  in  Fig.  6.  This 

block  diagram  is  "equivalent"  in  the  sense  that  its  reduction  via  loop  closures 

and  block  diagram  algebra  results  in  Eq  49  or  50  for  the  open  command  (q  )  loop, 

and  in  Eq  kj  or  48  for  the  complete  closed-loop  system.  In  one  limiting  case, 

when  G12  =  0  (the  q2  loop  opened),  the  resulting  single-loop  system  is  just  the 

command  loop  closed  around  the  vehicle-alone  transfer  function,  Q1c  .  For  the 

1 02 

other  limiting  case,  G21  =  0,  the  q2  loop  is  the  only  one  closed,  and  the  modi¬ 
fied  vehicle  dynamics  are  those  resulting  from  the  closure  of  a  single  unity- 
feedback  loop  about  the  transfer  function  G12Q2^  . 

A  very  instructive  way  to  consider  the  effect  of  the  closed  loop  is  to 
lump  all  of  its  consequences  into  "effective"  vehicle  changes.  In  this  view, 
closure  of  the  q^  loop  changes  the  command  loop  effective-vehicle  transfer 
function  from  Q1  to  (see  Fig.  6).  Referring  to  Eq  49  and  50,  the 

effective-vehicle  characteristics,  Q-|  will  be  the  open- loop  transfer  function 
of  the  command  loop  with  the  controller  transfer  function  G21  removed. 
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In  general. 
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>ck  Diagram  for  the  System 


Therefore  the  so-called  root  locus,  or  high  frequency,  gain  of  Q] is  identical 


to  that  for  Qi 


52- 


The  primed  notation  on  Qi' gg  is  used  to  indicate  merely  that  one  loop  has 

been  closed,  and  is  not  intended  to  specify  the  particular  loop  closure  involved. 

If,  for  example,  a  loop  closure  of  q^— fc-S-j  replaced  the  qg— *-S-|  closure  in 

Fig.  5,  the  transfer  function  Q.'  would  still  indicate  the  resulting  command 

1  o2 

loop  effective-vehicle  transfer  function.  The  precise  meaning  of  the  primed 
notation  consequently  depends  on  the  local  context.  (Later,  the  primed  notation 
is  used  in  a  similar  fashion  on  individual  transfer  function  terms  to  indicate 
the  number  of  prior  loop  closures.) 

Finding  the  complete  form  of  Q-J g  involves  operations  converting 
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Such  operations  are  easily  accomplished  by  closing  the  single  loops  shown  in 
Fig.  7  using  any  of  a  'variety  of  conventional  servoanalysis  techniques.  Letting 
Gi2  =  N12/D12,  where  both  N-j  g  and  D-j  2  are  polynomials  in  s,  with  the  leading 
coefficient  of  D-j  g  equal  to  unity  and  the  leading  coefficient  of  K-j2  being  the 
root  locus  gain  of  G-jg,  the  results  desired  from  the  loop  closures  are: 
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(Closed-loop  pole  factors  of  (T) ) 
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a)  Closure  (T) 


b)  Closure  (2) ,  "Coupling  Loop" 


Figure  7.  q2  Loop  Closures  Involved  in  the  System 
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where  ©  and  ©  refer  to  the  single-loop  systems  shown  in  Fig.  7*  For  fixed 
vehicle  characteristics,  G-|  2  is  the  only  variable  quantity  in  both  relationships. 
Thus,  choosing  G-|  2  appropriate  to  either  ©  or  automatically  determines  the 
other  closure  (©  or  ©)  .  The  closures  of  ©  and  ©  are,  therefore, 
"simultaneous."  Loop  ©  is  a  primary  closure  directly  affecting  the  vehicle's 
poles,  whereas  loop  ©  (as  will  be  seen  below)  affects  the  vehicle  zeros  perti¬ 
nent  to  a  particular  input  control  deflection — in  this  case,  62.  Accordingly, 
loop  <2>,  or  its  reflection  in  Fig.  6,  is  often  called  a  "coupling  loop." 

The  two  loops  are  further  related  in  that  the  high  frequency  open-  and  closed- 
loop  asymptotes  of  systems  (T)  and  (2)  are  usually  nearly  identical,  i.e.,  the 
"root-locus  gains"  are  almost  the  same.  This  is  easiest  to  see  for  typical  air¬ 
craft  characteristics  in  which  the  highest  order  s  terms  in  A  stem  from  the  main 
diagonal  (an  822*33)  >  and  the  control  effectiveness  terms  (Fij)  all  have  the  same 
form  in  s  (usually  constants).  Then,  at  frequencies  much  greater  than  all  the 
poles  and  zeros  in  A,  denoted  as  "|s|  large",  pertinent  quantities  become 
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and,  similarly. 
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When  Fq-]  »  fiif22/f12?  a  relatively  common  occurrence  in  aircraft,  Eq  59  and  60 
are  very  nearly  the  same.  This  correspondence  should  be  kept  in  mind  when 
actually  performing  the  simultaneous  closures  of  ©  and  . 

When  the  results  of  the  simultaneous  closures  are  combined  and  multiplied 
by  the  vehicle  transfer  function  Qig^,  then  Qi'gg,  the  effective  q-| /&2  transfer 
function  of  the  vehicle  with  the  qg  loop  closed,  is 
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S<Flg2  (Closed-loop  pole  factors  of  (2)) 
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As  already  noted  in  connection  with  Eq  55  and  56,  the  net  "root-locus  gain"  for 
Q.|  ,  after  the  closures  of  (7)  and  (2),  will  be  the  same  as  that  for  Qi&0* 

This  high-frequency  gain  is  implicitly  contained  in  the  Q-f numerator  term 
resulting  from  closure  of  loop  (2)  (it  arises  from  the  Di2Nql52  numerator  term 
in  Eq  58) • 


After  Qj52 
are  determined 


is  found,  the  final  q-j/q^  closed-loop  system  characteristics 
by  another  loop  closure  using  the  system  of  Fig.  8. 


55 


Figure  8.  Command  Loop  Closure  for  "the  System 
q1c  &2  j  12  — B1 


In  summary,  the  steps  Involved  in  the  analysis  of  the  system  in  Fig.  5  are 

1 .  The  control  channels  are  divided  into  two  categories,  inner 
and  "outer"  loops,  reflecting  the  closure  sequence.  The  q2 
(or  5i )  loop  was  closed  first,  and  thus  was  the  "inner"  loop; 
whereas  the  q^  (or  52)  loop,  being  closed  second,  was  the 
"outer"  loop. 

2.  The  "inner  loop,"  Gi2Q2Qi  ((D>  7)  is  closed  with 

tentatively  selected  equalization  and  gains,  and  the 
closed  inner-loop  roots  are  found.  These  roots  become 
the  vehicle's  poles  for  the  outer- loop  closure. 

3.  Using  the  same  gain  and  equalization  selected  above, 
i.e.,  the  same  G-]2>  the  coupling  loop  (©,  Fig.  7)  is 
closed.  The  closed-loop  roots  resulting  from  this 
closure  become  the  vehicle's  zeros  for  the  outer- loop 
closure . 

k.  The  outer  loop  is  closed  in  a  conventional  manner  around 
the  modified  outer-loop  vehicle  transfer  function. 

5.  Possible  repetitions  of  steps  2  through  4  with  different 
equalizations  and  gains  if  the  result  of  step  4  is  not 
satisfactory. 

The  generalization  of  these  steps  to  handle  more  complex  systems,  such  as  that 
of  Fig.  4,  is  a  straightforward  extension. 
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E.  APPLICATION  TO  SYSTEM  SYNTHESIS 


Fundamentally,  the  five  analysis  steps  summarized  above  can  represent  one 
cut  in  a  cut-and-try  synthesis  procedure.  In  such  procedures  a  trial  system  is 
analyzed  to  determine  its  dynamic  characteristics,  which  are  then  compared  with 
dynamic  performance  objectives.  Deficiencies  revealed  by  the  comparison  are, 
hopefully,  eliminated  or  reduced  by  modifications  resulting  in  a  new  system  for 
trial.  This  system,  in  turn,  is  analyzed  and  assessed,  etc.  The  iterations 
continue  until  the  trial  system  characteristics  are  consonant  with  the  dynamic 
performance  objectives. 

The  number  of  iterations  required  to  achieve  the  penultimate  system  depends 
on  the  designer's  intuition,  on  his  ability  to  transfer  insights  into  analytical 
procedure,  and  on  his  capacity  to  draw  new  knowledge  and  understanding  from  the 
analysis  results.  An  analysis  procedure  cannot  substitute  for  background  and 
experience — -but  the  procedure  can  be  matched  with  a  presumed  background  so  as 
to  achieve  a  balanced  analyst-analysis  procedure  "system."  The  background 
presumed  as  a  match  for  the  multiloop  analysis  method  advanced  here  is:  an 
intimate  knowledge  of  single-loop  servoanalysis  techniques;  a  detailed  under¬ 
standing  of  vehicle  dynamics;  and  a  thorough  appreciation  of  the  changes  in 
effective -vehicle  dynamics  caused  by  idealized  (i.e.,  no  sensor  or  servo  lags) 
single-loop  controllers.  Assuming  this  core  technology,  the  multiloop  analysis 
procedure  has  been  designed  to  maximize  the  insightful  generating  aspects  and 
clarity  of  each  analysis  step. 

When  full  advantage  is  taken  of  the  insights  provided  by  the  presumed 
background  level,  the  multiloop  analysis  procedure  can  form  the  basis  for 
almost  direct  synthesis  (i.e.,  the  number  of  analysis  iteration  cycles  required 
to  achieve  the  final  system  approaches  one).  The  trick,  of  course,  is  to  start 
off  with  a  good  trial  system,  which  amounts  to  almost  knowing  the  answer.  To 
illustrate  the  possible  impact  of  background  information  on  the  evolution  of 
the  initial  trial  system  a  few  items  are  listed  below,  using  the  two-loop 
system  of  Fig.  as  a  particular  example. 

1 .  Hints  about  desirable  inner-loop  (qg  loop)  characteristics 
can  be  revealed  by  studies  of  the  outer- loop  system  using 
several  possible  alternatives  for  the  vehicle-alone  trans¬ 
fer  function,  Q-j  g  ,  as  representative  of  Q-Jg  .  In  such 
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studies,  the  approximate  effects  of  possible  inner-loop 
feedbacks  on  Qi  g2  can  be  determined  readily  using  the 
literal  approximate  factors  of  Qi  g2  (see  Refs  2  and  3) 
and  the  equivalent  stability  derivatives  corresponding 
to  idealized  controllers  (Refs  2—5) . 

2.  An  appreciation  for  the  approximate  forms  obtainable  for 
the  closed  inner  loops  and,  consequently,  for  the  outer- 
loop  effective-vehicle  transfer  function  Q-fgg,  can  136 
advanced  from  studies  of  the  inner- loop  closure  involving 
Q2q  .  Again  approximate  factors  and  equivalent  stability 
derivatives  as  well  as  single-loop  servoanalysis  are  used. 

3.  Preliminary  indications  of  tradeoffs  between  parallel 
(the  q2  loop  closure)  and  series  (G21 )  equalization  are 
implicit  in  the  above  single-loop-closure  studies.  Thus, 
single-loop  closures  about  the  various  Qi  g2  forms  con¬ 
sidered  indicate  the  corresponding  types  of  equalization 
needed  in  G21  ;  and  the  approximate  forms  possible  (with 
particular  inner- loop  feedback  quantities)  and  desirable 
(with  particular  forms  of  G21  )  for  Ql' g2  can  be  inferred. 

Studies  of  the  above  types  either  involve  vehicle  characteristics  as 
variable  or  emphasize  controller  equalization  as  variable.  Both  types  of 
studies  use  well-known  analytical  operations  and,  for  aircraft,  a  large  body  of 
quite  general  results  exist  (see,  for  example.  Refs  2,  3,  1 3,  and  1^).  Using 
the  approximate,  but  highly  indicative,  information  gathered  from  such  past 
results  as  a  guide,  and  adding  the  additional  analytical  complication  required 
to  account  for  the  dynamics  of  such  elements  as  sensors  and  servos,  an  excellent 
first-cut  trial  system  is  easily  obtained.  The  exact  analysis  routine  summa¬ 
rized  previously  can  then  be  used  in  the  final  steps  of  a  practical  synthesis 
procedure.  For  systems  as  simple  as  that  shown  in  Fig.  5>  repetitious  opera¬ 
tions  can  often  be  avoided  completely  by  this  bringing  of  previous  single— loop 
system  knowledge  to  bear  upon  the  multiloop  problem. 

In  spite  of  the  above  remarks  cn  synthesis,  it  is  well  to  remember  that  the 
process  outlined  does  not  involve  unique  operations  in  compliance  with  a 
straightforward  routine.  One  of  the  more  tricky  aspects  has  been  casually 
bypassed  until  now.  This  is  the  selection  of  a  particular  block  diagram,  or 
sequence  of  closures,  from  the  several  possible.  The  block  diagram  itself  has 
little  significance — it  is  the  loop  closure  sequence  represented  graphically  by 
the  block  diagram  which  is  important  in  eliminating  or  reducing  iterations  when 
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the  analysis  technique  is  use!  as  a  synthesis  tool.  From  the  point  of  view  of 
pure  analysis  the  closure  sequence  is  immaterial.  But  in  synthesis  the  closure 
sequence  can  be  all-important.  For  instance,  certain  loops  are  necessarily 
closed  before  others  for  which  they  provide  parallel  equalization;  and  the 
necessary  use  of  incomplete  loop  closure  criteria  causes  some  loop  closure 
sequences  to  result  in  extreme  variations  in  loop  adjustments  as  a  function  of 
iteration,  whereas  for  other  sequences  there  is  practically  no  change  from  one 
iteration  to  the  next,  etc.  Because  of  considerations  like  these  there  are  at 
least  some  preferred  loop  closure  sequences  and  sometimes  even  a  uniquely  desired 
one.  Desirable  sequences  are  not  always  simple  to  determine.  Fortunately,  for 
vehicular  control  systems  a  set  of  factors  can  be  promulgated  which,  when 
properly  considered,  will  ordinarily  provide  the  insight  needed  to  construct  a 
unique  block  diagram  (or  closure  sequence)  which,  in  the  practical  sense  of 
minimizing  interloop  interactions,  will  also  minimize  iterations.  These  are 
summarized  in  Table  I. 

As  should  be  apparent  from  Table  I,  many  of  the  detailed  insights  involved 
can  stem  from  analysis  of  single-loop  systems  and  from  sinqplifled  multiloop 
analyses  -  Thus,  the  equivalent  stability  derivative  approach,  vehicle  approx¬ 
imate  transfer  functions,  and  single -sensor-loop  studies  already  mentioned 
play  yet  another  dominant  role  in  multiloop  vehicular  systems  synthesis. 

After  thorough  consideration  of  the  factors  in  Table  I,  a  quasi-unique 
block  diagram,  or  sequence  of  closures,  can  usually  be  established.  The  next 
problem  is  that  of  closure  criteria  for  .the  several  loops  .  Again  only 
engineering  judgment  factors,  which  must  ultimately  be  translated  into  concrete 
performance  measures,  can  be  delineated.  The  factors  to  be  considered  in 
establishing  the  actual  closure  criteria  for  each  loop  include: 

1  .  The  use  of  an  inner  loop  as  equalization  for  a  subsequent 
outer  loop . 

2.  Stability  and  response  of  loops  which  may  be  outer  loops 
in  one  mode  of  operation  and  inner  loops  in  another 
operational  mode . 

3-  The  "sensitivity  factors"  (pole,  zero,  and  gain  sensi¬ 
tivities,  Ref  10)  of  the  various  loops  and  their 
influence  on  the  outer  loops  .  Ordinarily,  inner-loop 
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closures  selected  should  he  such  that  the  effects  of 
inner- loop  parameter  variations  on  the  outer  loops  are 
as  small  as  possible. 

These  considerations,  tales  in  content  with  over-all  system  performance  specifi¬ 
cations,  can  ordinarily  serve  as  guidelines  for  the  selection  of  specific  closure 

criteria. 
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SECTION  III 


LONGITUDINAL  EXAMPLE  -  ALTITUDE  CONTROL  SYSTEM 

To  illustrate  the  previous  general  treatment  two  concrete  examples  will  he 
presented,  one  for  aircraft  longitudinal  control  (in  this  section)  and  the 
other  for  lateral  control  (in  Section  IV)  .  Loop  closure  sequence  and  criteria 
constitute  a  large  part  of  the  practical  synthesis  problem,  so  their  selection 
wj_]_X  receive  detailed  attention  in  the  ensuing  treatment.  In  particular, 
the  general  factors  involved  in  the  selection  of  loop  closure  sequence  (Table  I) 
will  be  specialized  for  each  of  the  examples;  and  some  of  the  background  con¬ 
siderations  pertinent  to  closure  criteria  will  be  discussed  prior  to  the  selec¬ 
tion  of  representative  quantities  and  values .  Once  the  sequence  and  criteria 
are  fixed  the  analysis  itself  is  routine,  and  the  detailed  examples  are  then 
presented  as  straightf orward  numerical  exercises . 

The  longitudinal  system  is  an  altitude  control  system  wherein  the  airframe 
is  controlled  by  the  elevator  which,  in  turn,  is  activated  by  feedbacks  involv¬ 
ing  the  pitch  angle,  0,  and  the  altitude,  h.  This  system  is  illustrated  in  the 
block  diagram  of  Fig.  9.  Here  the  usual  "e"  subscript  on  the  elevator  deflec¬ 
tion  symbol,  6e,  is  eliminated  to  simplify  the  notation. 


Figure  9.  Longitudinal  Control  System 


This  system  has  a  very  simple  multiloop  form  because  only  one  control 
deflection,  5,  is  used.  The  surface  actuator  is,  therefore,  common  to  both 
feedback  channels — a  feature  emphasized  in  the  block  diagram  by  indicating  an 
actuator  transfer  function  block,  Gg,  separate  from  the  altitude  and  pitch 
channel  equalizations. 

A.  LOOP  CLOSURE  SEQUENCE  AND  GENERAL  CLOSURE  CONSIDERATIONS 

For  this  system  a  unique  sequence  of  closures  is  easy  to  justify  using 
Table  II,  as  a  special  case  of  Table  I,  to  illustrate  the  reasoning.  Some  of 
the  less  obvious  remarks  in  Table  II  may  become  more  apparent  after  the  discus¬ 
sion  below.  All  of  the  factors  considered  in  Table  II  indicate,  or  are  compat¬ 
ible  with,  a  closure  sequence  with  h  as  the  outer  loop.  The  block  diagram  of 
Fig.  9  can  now  be  shown  with  an  hc  command,  as  in  Fig.  10,  and  the  system  can 
be  designated  as 

hc,  0— —5  (62) 


Figure  1 0.  Altitude  Control  System 
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B.  SPECIALIZATION  OF  SYSTEM  EQUATIONS 


The  controller  equation  for  the  system  is 


S  -  GBh(hc  -  h)  -  G500 


(63) 


where  Gg^  -  G^Gg  and  GBq  -  GgGg.  As  a  special  case  of  the  general  equations 
in  the  last  section  the  quantities  involved  here  are  as  given  below. 


Control  Deflections : 


51  =  0 

52  =  5 


Output  -  Motion  Quantities : 


Vehicle  Transfer  Functions : 


Controller  Transfer  Functions: 
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The  closed-loop  transfer  function  h/hc  can  be  determined  directly  from  Eq  39 
and  40  if  Fig.  10  is  considered  a  special  case  of  Fig.  4  and  the  terms  and 
symbols  adjusted  accordingly. 


jTj  _ G8hNhs _ 

hcJh,e-—  5  A  +  G50%5  +  G8hNhB 

_  G5hHS _ 

1  +  G8^S  +  G8hHB 


(64) 


For  the  purposes  of  this  report  all  literal  airframe  transfer  functions  are 
written  as  in  Refs  2,  3,  or  11,  all  of  which  are  compatible.  Numerical  trans¬ 
fer  functions,  when  required  for  plotting  or  other  purposes,  are  based  on  the 
typical  case  given  in  Ref  11 .  However,  although  derived  using  these  numerical 
values,  many  plots  are  presented  in  generic  form  with  the  poles  and  zeros 
identified  by  their  literal  values  to  enhance  the  clarity  of  presentation.  In 
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terms  of  such  literal  values  the  generic  airframe -a lone  transfer  functions, 
expressed  in  the  root-locus  form,  are: 


%s  Ag(s  +  1/Tfll)(s  +  1/Teg) _ 

(s2  +  2£j/UpS  +  o|)(s2  +  2CSjPsps  +  cii,2) 

%5  Ah(s  +  I/T^Hb  +  1/Thg)(s  +  1/Th^) 

s(s2  +  2^ptnpS  +  o^)(s2  +  2£spu>sps+  cug|) 


(65) 


(66) 


As  already  noted,  approximate  factors  relating  the  poles  end  zeros  in  Eqs  65 
and  66  with  the  aerodynamic  and  inertial  parameters  of  the  airframe  are  given 
in  Refs  2  and  3 • 

The  appropriate  equivalent  block  diagram  for  Pig.  10  is  shown  in  Pig.  11  . 
Hg,  the  equivalent  vehicle  transfer  function  with  the  9  loop  closed,  is  also 
shown  in  this  figure. 


hi  -M'-  H§ 
8J0-8”  8  I  +  g80®8 


- - 


Figure  11.  Equivalent  Block  Diagram  for  Altitude  Control  System 

h  ,  9  — -  5 
c 
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C.  INKER  LOOP  CONSIDERATIONS 


1  .  Inner  Loop  Alone 


Establishing  the  closure  criteria  for  the  two  loops  is  somewhat  more 
involved  than  the  closure  sequence  selection.  Considering  the  fundamental 
vehicle  dynamics  and  possible  disturbances,  and  recognizing  that  the  attitude 
loop  will  probably  be  the  outer  loop  (i.e.  attitude  commands)  in  some  oper¬ 
ational  mode  of  the  flight  control  system,  the  9  loop  alone  should  be  capable  of 

( 1 )  Providing  good  attitude  regulation  of  the  vehicle  in 
the  presence  of  external  disturbances  (e.g.  gusts), 
and  both  short  and  long  term  vehicle  asymmetries  such 
as  sudden  changes  in:  c.g.  with  stores  release,  thrust 
moment  with  power  changes,  external  configuration  with 
actuation  of  auxiliary  surfaces  such  as  flaps  or  brakes , 
etc.;  and  slow  changes  in:  c.g.  with  fuel  consumption, 
flight  operating  conditions,  etc. 

(2)  Accomplishing  commanded  changes  in  attitude  with  good 
response . 

For  nominal  airframe  characteristics  these  aims  can  be  accomplished  using  a 
controller  with  lead  equalization,  and  a  crossover  frequency,  co0,  somewhat 
greater  than  the  short  period  undamped  natural  frequency.  Typical  characteris¬ 
tics  for  such  a  system  are  illustrated  in  the  sketches  of  Fig.  12.  The  con¬ 
ventions  established  in  Fig.  12,  for  identifying  the  meaning  of  the  various 
line -values  used  on  the  Bode  plots,  are  followed  in  all  remaining  plots. 
Neglecting  the  controller  lag  dynamics,  the  open  loop  transfer  function  in 
Bode  form  will  be 


G(s)  = 
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+  1 


(6?) 


where  Tg  is  the  lead  equalization  time  constant.  ('This  equalization  can  be 
developed  directly  from  a  9  sensor  using  passive  elements,  or  from  a  rate 
sensor  plus  a  8  sensor.)  The  actual  values  for  the  closed-loop  transfer 
function  factors  are  readily  found  by  the  decomposition  technique  using  the 
G(  ju>)  and  G(  -a)  Bode  plots  of  Fig.  12,  or  from  a  root  locus.  These  techniques. 


(c)  Root  Locus 

Figure  12.  Open-  and  Closed-Loop  Characteristics  of  Attitude  Control  System 
with  Ideal  Controller  and  Nominal  Airframe 


and  others  Involved  in  unified  servoanalysis  procedures,  are  covered  elsewhere 
(Ref  9). 

When  closed  with  values  of  open-loop  dc  gain,  K,  as  shown  in  Fig.  12, 
the  resulting  closed-loop  transfer  function  is 
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Eg  68  uses  several  aspects  of  the  notational  scheme  adopted  in 
this  report  for  multiloop  analysis  results,  so  a  short  aside  is 
pertinent  to  discuss  notation  using  this  concrete  example.  All 
of  the  closed-loop  poles  have  one  prime,  indicating  that  they 
result  from  one  loop  closure.  Also  the  actual  notation  selected 
for  the  closed-loop  poles  reflects  their  approximate  values .  Thus 
l/T q' ,  l/T q'  and  l/T£,  are  a11  near  the  open-loop  numerator  terms 
1/Tq] ,  l/T q  ,  and  l/TE,  respectively.  If  a  much  lower  value  of  K 
had  been  used,  such  that  the  closed-loop  factors  remained  quadratics 
(i.e.  K  <  K-| ) ,  then  the  closed-loop  denominator  would  have  had  the 
form 
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In  this  case  the  notation  used  for  the  closed-loop  quantities 
reflects  their  origin  in  open-loop  quantities.  If  intermediate 
values  of  gain  had  been  used,  such  that  the  closed-loop  phugoid 
were  more  than  critically  damped  but  the  closed-loop  short  period 
were  still  quadratic  (i.e.  K-,  <  K  <  K2),  then  the  closed- loop 
denominator  would  have  the  form 
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So,  as  gain  is  increased,  the  notation  for  the  closed-loop 
phugoid  and  short  period  modes  has  the  following  transitions : 
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Phugoid : 

Open-loop  (K  =  0) 


Closed-loop,  quadratic  (K  <  K-|  ) 

Closed-loop,  overdanrped  (K-j  <  K) 
Short  Period: 

Open-loop  (K  =  0) 


Closed-loop,  quadratic  (K  <  iC,) 

Closed-loop,  overdamped  (K^  <  K) 

The  rules  for  this  notational  sequence  can  he  summarized  as  : 

1  .  The  number  of  primes  present  indicates  the  number 
of  loop  closures . 

2 .  When  the  closed-loop  form  is  the  same  as  the  open- 
loop  form,  the  notation  for  the  closed-loop  factor 
is  the  same  as  the  open-loop  factor,  plus  a  prime. 

In  this  case  the  origin  of  the  closed-loop  factor 
is  always  at  hand. 

J.  The  notation  for  closed-loop  factors  which  differ 
in  form  from  their  open-loop  origins  reflects  the 
open-loop  zeros  which  the  factors  approach  as 
K — •—  When  no  open-loop  zero  exists  in  this 

situation,  a  special  form  is  coined  which  reflects 
the  origin  of  the  fact  (e.g.  l/Tgp). 

These  rules  are  special  cases  of  those  summarized  in  the  nomen¬ 
clature  section  at  the  front  of  this  report. 

For  the  attitude  system  closed  with  high  gain,  as  in  Fig.  12,  the 
closed-loop  dc  gain,  K/(l  +  K) ,  is  approximately  one,  and  the  dipole  pairs 
(Tq1  s  +  l)/(T^s  +1),  (Tq2s  +  D/(Tg2s  +  1),  and  (TEs  +  1)/(t£s  +  l)  in  e/ec 
all  nearly  cancel.  Thus  Eq  68  becomes  approximately 
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As  an  attitude  command  system,  this  is  very  good  indeed;  also  attitude  regulation 
properties  are  adequate.  So  from  the  standpoint  of  attitude  control,  the  system 
shown  is  satisfactory.  However,  two  caveats  axe  worth  mentioning.  First,  for 
some  flight  conditions  the  low-frequency  region  centered  about  oip  may  not 
exhibit  such  large  amplitude  ratio  relative  to  that  near  short  period,  thereby 
making  it  difficult  to  attain  values  of  K  much  greater  than  one.  Such  condi¬ 
tions  can  be  corrected  by  low  frequency  lag-lead  equalization  (e.g.  adding  an 
integral  of  9e  to  the  controller).  Second,  the  controller  assumed  is  ideal, 
with  servo  and/or  sensor  lags  ignored,  and  the  loop  gain  can,  theoretically,  be 
increased  indefinitely  without  any  stability  difficulties .  Controller  lags  will 
be  taken  into  account  later  in  the  detailed  analysis.  Their  presence  limits  the 
allowable  gain  for  stability,  and  is  one  reason  lead  equalization  (i.e.  the 
f  s  +  1  factor)  is  needed  in  the  attitude  loop. 

JIj 

2.  The  Inner  Loop  as  Equalization  for  the  Outer  Loop 


When  considered  in  company  with  the  h  outer  loop  the  attitude  loop 
closure  should  be  of  such  form  as  to  relieve  outer  loop  equalization  require 
ments .  With  the  ideal  controller  considered  above,  the  effective  airframe 
altitude  transfer  function  Hg  becomes 
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The  forms  of  H&  and  H£  are  illustrated  in  the  sketches  of  Fig.  1  J.  To 
emphasize  the  basic  differences  in  the  dynamic  rather  than  static  forms,  a  common 
gain  factor  is  used.  ""Actually  H5  would  appear  (l  +  KgJCgg)^  =  (KgKeg)^  below 
the  plot  for  (1  +  KgKgg)Hg.  Comparison  of  these  forms  (which  are  minimum  phase 
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Figure  1  3.  Comparison  of  Hg  and  (l  +  KeKeB)Hg, 


except  for  the  breakpoint  at  1 /T hj)  readily  indicates  the  superiority  of  H5  over 
H5.  For  a  pure  gain  altitude  controller  operating  on  Hg,  the  crossover  frequency 
must  be  well  below  l/Tj^,  (to  avoid  the  lightly  damped  phugoid)  and  the  resulting 
closed-loop  system  response  will  exhibit  a  very  lightly  damped  or  even  slightly 
unstable  oscillation  of  the  phugoid  variety.  When  the  vehicle-alone  phugoid  is 
unstable,  as  it  may  be  for  flight  conditions  involving  pitching  moment  changes 
with  speed,  a  stable  pure  gain  closure  is  impossible.  On  the  other  hand,  a 
pure  gain  h  controller  is  completely  suitable  for  the  system,  and  stable 
crossovers  in  the  neighborhood  of  1 /Tg  are  readily  obtained.  The  dominant 
closed-loop  mode  is  then  a  quadratic  having  an  undamped  natural  frequency  near  cq,, 
and  the  gain  can  be  selected  to  yield  well-damped  responses.  For  these  statements 
to  be  true ,  the  inner  loop  gains  must  be  high  enough  to  overdamp  the  phugoid, 
and  preferably  high  enough  to  force  the  l/Tg  pole  to  the  neighborhood  of 
l/T a2,  thereby  effectively  maximizing  available  bandwidth.  High  gain  Q  inner 
loops  are,  therefore,  a  fundamental  requirement  for  pure  gain  altitude  controllers. 

An  obvious  alternate  to  the  above  system  is  to  use  lead  equalization 
in  the  h  loop.  This  is  sometimes  necessary  anyway  (e.g.  in  terrain  following), 
hut  it  always  introduces  additional  difficulties.  For  example,  the  system 
(•with  no  inner  attitude  loop)  will  produce  drastic '  pitch  angle  changes  when 
the  airframe  hits  a  vertical  gust.  Considerations  of  this  nature  are  important 
when  comparing  competing  systems  for  the  altitude  control  task. 

A  final  point  worth  mentioning  is  the  possibility  of  1  /Thl  becoming 
negative .  As  discussed  in  Ref  2  this  depends  primarily  upon  the  drag  and 
thrust  characteristics  of  the  vehicle.  In  practice  negative  values  of  1 /Thl 
cannot  he  avoided  near  the  absolute  ceiling  or  in  landing-approach  when  flying 
at  speeds  below  that  for  minimum  drag.  When  this  occurs  an  altitude  system 
using  only  elevator  control  will  exhibit  a  divergence  (with  a  time  constant 
which  approaches  Thi  )  and  will  therefore  be  unacceptable. 

D.  NUMERICAL  EXAMPLE 

The  preceding  discussion  covers  most  of  the  general  aspects  of  an  altitude 
control  system  with  an  inner  attitude  loop.  The  present  article  provides  a 
numerical  example . 
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To  make  the  example  more  realistic,  a  high  frequency  second-order  lag  is 
used  to  approximate  the  servo  characteristics,  and  a  first-order  lag  to  simu¬ 
late  altitude  sensor  dynamics.  Thus,  the  servo  transfer  function  is 
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The  inner-loop  controller  transfer  function  is. 
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and  the  outer-loop  controller  transfer  function  becomes 
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where  Ta  is  the  time  constant  of  the  altimeter  installation. 
Numerical  values  for  the  altitude  sensor  and  servo  will  be  taken  as 

5b  =  0-7,  ^  =  50, 
l/Ta  =  15 


(74) 


The  vehicle  characteristics  are  similar  to  the  nominal  case  previously 
used  as  an  example  throughout  Ref  1 1  .  These  are : 
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The  attitude  channel  equalization  time  constant  would  ordinarily  be  found 
by  an  attitude  subsystem  optimization  procedure  which  considered  the  vehicle 
dynamics  over  the  entire  flight  regime.  This  is  beyond  the  present  scope,  so 
a  value  of  l/Tg  =  2.4,  which  is  reasonably  representative,  will  be  arbitrarily 
selected  for  the  example. 

With  all  required  numerical  values  fixed  the  analyses  can  proceed  in  a 
straightforward  fashion  to  determine  the  two  gains  KQ  and  K,  .  The  first 
objective  is  to  find  the  poles  and  zeros  of  l/(l  ■+■  Gg0©^)  .  Siggy  and  job  Bode 
plots  of  the  open-loop  transfer  function  Gg0©g  are  presented  in  Fig.  14,  and 
a  root  locus  sketch  (not  to  scale)  is  shown  in  Fig.  13.  The  actual  closure 
used  depends  in  practice  on  several  factors  (e.g.  whether  an  adaptive  device 
is  to  be  used,  effects  of  parasitic  nonlinearities,  etc.)  over,  and  above  those 
general  considerations  discussed  in  the  second  article  of  this  section  (i.e. 
bandwidth,  response  of  9  loop  to  commands,  etc.).  To  make  matters  simple,  a 
closure  criterion  based  on  a  phase  margin  of  40  degrees  is  used  in  Fig.  14. 

The  resulting  system  is  compatible  with  all  the  general  considerations,  i.e. 
it  will  exhibit  good  transient  response  to  attitude  commands,  good  attitude 
regulation,  reasonable  insensitivity  to  likely  parasitic  nonlinearities,  etc. 
However,  the  criterion  may  not  be  appropriate  for  certain  types  of  adaptive 
gain-changing  schemes,  especially  those  involving  low  amplitude  limit  cycles 
as  an  essential  characteristic . 

For  a  phase  margin  of  40  degrees  the  gain  K0  will  be  8.3  db  or  2.6  in 
linear  units  .  Using  the  decomposition  and  other  techniques  of  the  unified 
servoanalysis  method  (Ref  9),  the  closed-loop  transfer  function  G&0@g/(l  +  Gg0©&) 
becomes  (See  Fig.  14  -  dotted  line) 
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Closed-Loop  Bodes  for  Example  Attitude  Control  Subsystem 


Figure  15.  Root  Locus  Sketch  for  Example  Attitude  Control  Subsystem,  Gse®5 


Eq  77  reveals  that  closing  the  loop  has  overdamped  the  phugoid  and  short  period 
modes  and  reduced  both  damping  ratio  and  undamped  natural  frequency  for  the 
quadratic  mode  associated  with  the  servo.  The  closed-loop  phugoid  factors  are 
quite  close  to  the  zeros,  1/T9i  and  l/T02;  and  one  of  the  closed-loop  short 
period  factors  (l/T^)  is  near  the  1  /Tg  zero.  Thus,  the  low  frequency  poles  and 
zeros  are  largely  self -canceling  and  the  approximate  closed-loop  transfer  func¬ 
tion  for  the  attitude  loop  is 
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Eq  78  is  the  equivalent  of  Eq  69,  with  the  presence  of  the  servo  dynamics 
giving  rise  to  the  additional  second  order  factor. 

The  closed-loop  attitude  system  transfer  function  of  direct  interest  for 
the  subsequent  altitude  loop  closure  is 
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This  transfer  function,  when  combined  with  the  airframe  transfer  function  Hg, 
gives  Hg,  the  effective  airframe  altitude  transfer  function  with  the  attitude 
loop  closed,  that  is. 


Except  for  the  quadratic  pair  due  to  the  servo,  Eq  80  is  identical  in  form  to 
Eq  70,  which  was  derived  for  an  ideal  no-lag  controller.  When  the  altitude 
controller  transfer  function  Ggjj  is  combined  with  HJ,  the  numerator  servo 
quadratic  will  exactly  cancel  the  servo  characteristics  in  the  denominator 
of  Ggjj.  The  effective  servo  dynamics  for  the  h  loop  closure  therefore  result 
from  the  attitude  closed-loop  system. 
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Adding  controller  dynamics  to  Eq  80  results  in  the  open-loop  altitude  system 
transfer  function. 
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The  jffi  and  siggy  Bode  plots  for  this  transfer  function  are  given  in  Fig.  1  6  and 
a  root  locus  sketch  is  shown  in  Fig.  17.  As  in  the  attitude  loop,  the  open- 
loop  characteristics  in  the  frequency  region  of  possible  crossovers  are  such 
as  to  result  in  good  closed-loop  dynamic  response  if  a  phase  margin  of  40  to  60 
degrees  or  so  is  used  for  the  closure  criterion.  For  a  phase  margin  of  40 
degrees,  Kh  will  be  -49.4  db  or  O.OO34  rad/foot. 

The  complete  closed-loop  system  poles  can  be  determined  readily  using 
unified  servoanalysis  techniques.  However,  the  entire  procedure  is  simplified 
considerably,  with  little  loss  in  accuracy,  if  advantage  is  taken  of  a  simplified 
form  of  G5hH£.  In  Fig.  16  it  will  be  noted  that  |GghH£|  is  much  less  than  one 
in  the  frequency  region  above  the  breakpoint  at  T/T£.  Thus,  the  contributions 
of  ^h11^  to  the  system  frequency  response  will  be  very  small  at  frequencies 
above  oj  =  l/T£  =3-5  and  these  contributions  will  be  changed  only  very  slightly 
if  the  amplitude  asymptote  through  the  break  point  at  CD  =  1  /Tq  =  1  .05  is 
taken  as  the  high-frequency  asymptote  of  Gg^.  Thus,  a  simplified  open-loop 
transfer  function 


(82) 


is  capable  of  exhibiting  the  major  changes  introduced  by  the  h  closure.  The 
high  frequency  portion  of  this  simplified  plot  is  shown  by  the  dashed  line  on 
Fig.  16.  The  higher  frequency  terms  not  contained  in  Eq  82  will  be  modified 
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only  slightly  from  those  given  in  Eq  8l  by  the  altitude  closure.  Using  the 
decomposition  technique  on  this  simplified  open-loop  transfer  function,  a 
simple  approximate on  to  the  closed-loop  transfer  function  is 
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Eq  83  gives  only  the  dominant  mode  for  the  altitude  control  system.  Including  an 
approximation  to  the  higher  frequency  terms  from  Eq  8l ,  the  total  system  closed- 
loop  transfer  function  would  he  approximately 
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E.  CONCLUDING  REMARKS 

With  the  control  gains  set  as  given  above  the  altitude  control  system  will 
exhibit  characteristics  generally  accepted  as  good.  The  transient  response  to 
an  altitude  command  will  be  rapid  and  well  damped,  and  the  shape  of  the  response 
to  a  step  or  cut-off  ramp  will  be  essentially  that  of  a  unit-numerator  second- 
order  system  having  a  damping  ratio  of  about  0.5  and  an  undamped  natural 
frequency  of  1  .0  rad/sec.  The  steady-state  altitude  error  will  be  zero; 
however,  the  airspeed  will,  in  general,  change  to  make  this  possible.  (To 
also  keep  the  airspeed  constant  would  require  a  more  sophisticated  system 
incorporating  some  other  control,  such  as  a  throttle  loop.)  The  effects  of 
vertical  gusts  will  be  largely  suppressed  by  the  attitude  inner  loop,  which 
will  also  tend  to  minimize  the  effects  of  both  sudden  and  slowly  varying  air¬ 
frame  changes  (e.g.  flap  deflections,  eg.  shifts,  etc.).  Further,  the  attitude 
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loop  as  a  separate  entity  will  also  exhibit  rapid,  well -damped  responses  (to 
commands)  which  look  like  those  of  a  low-order  system.  Only  a  slight  dynamic 
droop,  (due  to  incomplete  cancellation  of  l/T-g  and  l/l^,  and  1  T and  1  T g^) , 
and  a  small  static  droop  or  error  (i.e.  S/Sc  =  0.93  instead  of  1 .0  as  s  -*  0) 
detract  from  almost  ideal  behavior  for  the  attitude-alone  system.  Thus,  the 
analysis  procedure,  as  illustrated  in  this  section,  has  resulted  in  a  rapid 
and  relatively  straightforward  nominal  system  synthesis  without  any  repetitious 
analys is . 

There  remains  some  question  as  to  the  detailed  selection  of  criteria  and 
equalization  made  in  the  course  of  the  analysis.  Such  questions  can  easily  be 
answered  by  making  a  sensitivity  analysis  on  the  final  closure  results.  For 
example,  the  sensitivities  of  the  closed-loop  poles  to  changes  in  attitude  and 
altitude  gains,  attitude  equalization  time  constant,  etc.  (computed  using  the 
methods  of  Ref  10)  can  be  plotted  as  vectors  on  a  closed-loop  root  diagram  to 
illustrate  the  first  order  effects  of  such  changes.  These  vectors  indicate  the 
direction  and  magnitude  of  shifts  in  the  closed-loop  poles  due  to  specified 
per-unit  shifts  in  the  various  open-loop  quantities .  Thus ,  the  effects  of 
small  changes  from  the  nominal  values  used  in  the  present  analysis  would  be 
immediately  available  from  the  sensitivity  results. 

A  similar  system  synthesis  could  be,  and  normally  is,  accomplished  using 
the  analog  computer  as  a  tool.  Somewhat  more  trial  and  error  is  involved,  and 
some  sensitivities  are  a  bit  more  difficult  to  obtain.  (The  calculation  of 
sensitivities  is  a  simple  additional  step  when  the  system  analysis  data  used 
in  this  section  is  available.)  The  analysis  procedure  presented  here  was  not 
developed  to  be  competitive  with  analog  results  (although  it  is  superior  in 
many  ways) .  Instead  one  of  the  major  reasons  for  its  development  was  to 
achieve  an  alternate  pencil-and-paper  analysis  technique.  The  statement 
that  similar  results  would  be  obtained  in  an  intelligent  synthesis  program 
using  the  analog  computer  simply  indicates  that  the  desired  alternate  analysis 
method  has  been  achieved. 

As  emphasized  throughout  the  report,  the  loop  closure  sequence  is  extremely 
important  if  results  such  as  those  obtained  above  are  to  be  readily  achieved. 
The  altitude  control  system  provides  an  excellent  example  in  this  regard . 
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Consider,  for  the  moment,  that  the  h  loop  had  been  closed  first,  with  the  0 
loop  as  the  outer  loop.  Then  the  gains  derived  above  would  have  been  a  highly 
unlikely  result  because  of  the  involved  reasoning  required  to  establish  closure 
criteria.  For  example,  the  h  loop  would  have  to  be  closed  to  give  an  unstable, 
or  nearly  unstable  inner  loop,  which  would  then  be  stabilized  by  the  outer 
attitude  loop.  With  the  closure  sequence  actually  selected  the  criteria  were 
based  upon  relatively  simple,  essentially  "standard",  considerations  derived 
from  an  intimate  knowledge  of  the  characteristics  of  low  order  systems .  Now 
that  the  system  has  been  explored  the  sequence  is  not  nearly  so  important, 
since  special  criteria  could  now  be  evolved  to  allow  the  synthesis  of  a  reason- 
nble  system  even  with  the  altitude  loop  closed  first . 
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SECTION  IV 


LATERAL  EXAMPLE  -  BANK  ANGLE  CONTROL  SYSTEM 

As  noted  in  the  introduction  to  Section  III,  the  longitudinal  and 
lateral  numerical  examples  are  intended  both  to  illustrate  the  general  analy¬ 
tical  treatment  developed  in  Section  II,  and  to  describe  further  some  of  the 
considerations  leading  to  the  selection  of  loop  closure  sequence  and  criteria. 

The  insight  necessary  for  such  selection  can  be  developed  by  preliminary 
analyses  which  study  single-loop  closures  for  each  of  the  several  loops  involved, 
and  loop  closures  based  on  ideal  controllers  and  approximate  literal  vehicle 
transfer  function  factors.  In  the  longitudinal  example  single-loop  closures 
(of  Hg  and  loops)  were  discussed  and  the  results  were  used  to  guide  certain 
aspects  of  the  multiloop  analyses  .  The  lateral  controller  treated  below  will 
provide  an  example  using  the  other  technique  mentioned,  i.e.  approximate 
factors  and  ideal  controllers  converted  to  equivalent  stability  derivatives . 
Otherwise,  the  discussion  will  proceed  in  much  the  same  way  as  that  adopted 
for  the  longitudinal  case . 

A.  SPECIALIZATION  OF  SYSTEM  EQUATIONS 

The  lateral  system  chosen  as  an  example  is  basically  an  attitude  control 
system  intended  to  maintain  zero  yawing  velocity  and  zero  bank  angle.  It  can  alsc 
be  used  to  turn  the  aircraft  by  introducing  bank  angle  commands.  In  this 
system  the  rudder  is  activated  by  feedbacks  involving  yawing  velocity,  r, 
and  the  aileron  is  activated  by  the  bank  angle  error,  cp£.  The  block  diagram 
is  given  in  Fig.  1 8  and  the  system  designation  is 


r 


6 

6 


r 

a 


The  controller  equations  are 

Br  =  -°5rrr  =  "Grr 
5a  =  GSaT^c  ”  T)  =  Gcp(Tc  '  T) 


(85) 


(86) 
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Figure  18.  Bank  Angle  Control  System 
'Pc  — ~  5a>  r  8r 

where  the  conventional  dual  subscript  notation  for  the  controllers  (which 
specifies  both  output  and  input  controller  quantities)  can  here  be  simplified 
to  a  single  subscript  without  confusion.  The  system  is  basically  a  special 
case  of  the  generalized  system  shown  in  Fig.  5 .  In  terms  of  equations  in 
Section  II  the  quantities  involved  here  are: 


Control  Deflections : 


Output  Motion  Quantities : 


Vehicle  Transfer  Functions : 


<ll  =  <P 
<3p  =  r 


^16,  _ 

Ql82  =  $&a 

^25l  =  %r> 

Q262  =  R5a 
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Controller  Transfer  Functions : 


G1 1  -  G1  3  =  G22  =  G23  =  0 
G1 2  =  GSrr  =  Gr 
g21  =  GSacp  3  Gcp 


As  special  cases  of  the  closed-loop  transfer  function  forms  given  in  Eqs  47 
and  48, 


JL 

<PC 


+  ) 


-Jr— Sr  A  +  +  G^Nq^  +  G^g^) 

<P-~&a 


N, 


cp  r 


5*5-, 


GCP^6jl  +  Gr 


a 


cp  r 

,  %a5r, 

1  +  Gr%r  +  (l  +  Gr  -t^) 


>W 


(87) 


(88) 


The  open-loop  transfer  function  for  an  outer  cp  loop,  (presuming  the  r  —  8r 
channel  is  used  as  an  inner  loop  and  is  already  closed)  is  a  special  case  of 
Eqs  49  and  50, 


A  +  GpNpg 


(89) 


,  %yr, 

(1  +  Gr  -f^) 


N( 


G  - 

<P  5a  1  +  GrR5 


(  90) 


The  airframe  transfer  functions  in  literal  terms  will  be  taken,  for  this 
example,  to  have  the  forms. 


Nr5r  Ar(s  +  1/Tri)(s2  +  2Z]*±cs  +  a)2) 

A  (s  +  1/Ts)(s  +  l/TR)(s2  +  2^do)ds  +  cof) 
Ncpga  Aq,(s2  +  2(^ays  +  cd^) 

A  (s  +  l/Ts)(s  +  l/TR)(s2  +  2^d;nas  +  cof) 


(  9l) 

(92) 
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The  coefficients  of  these  transfer  functions,  and  literal  approximate  factors, 
in  terms  of  aerodynamic  and  inertial  parameters  are  presented  in  Refs  2  and  J. 
These  stem  from  the  three  degrees  of  freedom  equations 


which  describe  the  lateral  motions  of  the  aircraft  when  only  slightly  perturbed 

from  straight,  wings  level,  and  horizontal  flight.  By  analogy  with  Eq  54  the 

"coupling  numerator",  gr  ,  is  then 

bl  r 


0  Ygr  (s  -  Yv) 

<8rr  -  Ie»  %  -I i 

"4  "8r  -  H 

-  <4a»8r  -  * 

~  ^cpr^s  +  ^  /T<pr)  =  Kqjr(  Tqjrs  +  l)  (q4) 

Other  coupling  numerators  for  common  multiloop  systems  are  tabulated  in 

Appendix  B.  The  final  line  of  Eq  94  also  introduces  the  last  type  of  special 

notation  for  multiloop  analysis.  Unfortunately  even  the  double  subscripts  on 

K<pr  ana  T<pr  can  be  ambiguous,  although  not  for  most  practical  systems.  If 

essential  ambiguities  ever  occur,  a  sub  and  superscript  combination,  as  in 
cp  r 

%aSrj  could  conceivably  be  used  (although  a  special,  albeit  inconsistent, 
notation  would  probably  be  clearer) . 


Tl^ 


.  '  W] 

- 
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B.  Cp 


— -  Sa  AND  r  Sr  AS  SINGLE  LOOPS 


Because  two  control  deflections,  6r  and  6a,  are  used,  analysis  of  the 
system  will,  as  noted  above,  involve  a  "coupling  numerator”  closure .  The 
system  is,  therefore,  somewhat  more  complex  than  the  altitude  controller 
example  of  Section  III .  Another  general  difference  is  that  both  the  r  -*■  &r 
and  cp  — S&  feedbacks  must  play  an  attitude  stabilization  and  regulation  role, 
whereas  only  the  0  loop  had  such  requirements  in  the  longitudinal  system. 
Consequently  the  bandwidths  of  the  r  6  and  cp  — ■ —  6a  loops,  considered 
individually,  should  be  about  the  same  in  magnitude  Also  commands  into 
either  or  both  loops  will  provide  a  turning  capability  so  the  use  of  only  bank 
angle  commands,  cpc,  indicated  in  Fig.  1  8  is  somewhat  premature.  With  only  such 
rudimentary  considerations,  it  is  not  obvious  a  priori  which  should  be  the  inner 
loop.  More  insight  can  be  gained  by  considering  both  as  possible  inner  loops. 
Therefore  individual  loop  closures  of  r  6r  and  cp  —  6a  will  be  investigated 
below  before  delineating  the  loop  closure  sequence. 


1  .  r  -**■  Sr  as  a  Single  Loop 


The  Rq^  transfer  function  form  given  by  Eq  91  has  a  numerator  which  is 
full  of  surprises  if  all  manner  of  vehicles  are  considered.  The  cubic  can 
factor  into  either  a  first  order  plus  a  quadratic  (as  indicated  by  Eq  91 )  or 
into  three  first  orders.  Further,  the  roots  of  the  cubic  can  take  on  almost 
all  possible  variations  in  sign.  Therefore  a  complete  summary  of  the  pos¬ 
sibilities  for  Nr&  requires  a  large  number  of  different  conditions  to  be  taken 
into  account.  Some,  but  not  all,  of  these  are  codified  in  terms  of  approximate 
factors  in  Refs  2  and  15.  For  most  manned  aircraft,  however,  the  number  of  pos¬ 
sibilities  is  cut  down  substantially,  with  the  factors  usually  being  of  the 
following  nature : 
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T 


ri 


T  > 
R 


0 


COf.  « 


Oij.,  U>d  >  0 

?r  ^  0 
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I  £rl  «  1 


The  sketches  of  Fig.  19  illustrate  typical  dynamic  characteristics 
achievable  with  a  pure  gain  r  — -  5r  loop  closure  for  a  typical  manned  air¬ 
craft.  The  R§  transfer  function  has  the  features  listed  above,  and  both  £r 
and  l/T  are  taken  to  be  negative.  Fig.  19  indicates  that  the  spiral  can  be 
stabilized  readily  ( 1 /T^  is  positive)  and  that  considerable  damping  can  be 
added  to  the  Dutch  roll  mode  (£d  >  £d)  without  running  into  difficulties  due 
to  the  high  gain  stability  limitation  inherent  in  the  negative  value  assumed 
for  £r. 

The  same  conclusions  about  Dutch  roll  damping  augmentation  and  stabili¬ 
zation  of  the  divergent  spiral  mode  can  be  obtained  using  literal  approximate 
factors  and  the  method  of  equivalent  stability  derivatives .  As  given  by 
Refs  2  and  3,  the  literal  approximate  factors  for  the  spiral  and  Dutch  roll 
damping  are 


For  an  ideal  (no  lag)  controller  the  rudder  deflection  will  be 

S  =  -Krr 
r  r 


(97) 


which  will  create  a  yawing  acceleration  equal  to  -N§rKrr.  The  effective  value 
of  the  stability  derivative  ih.  will  then  become 


^effective  Nrvehicle  ^r^S-,. 

alone 


(93) 


The  effective  R^.  can  be  made  more  negative  (much  larger  in  magnitude)  than  the 
nominally  negative  R^  of  the  vehicle -alone.  Using  Rreff  to  replace  Rr  in  Eq  96, 
the  Dutch  roll  damping  (£d<%)  with  control  is  seen  to  be  increased  by  the 
increment  KrN&r/2.  Similarly,  the  substitution  of  Rreff  r°r  Nr  in  95  will 
increase  l/Ts.  For  values  of  Kr  such  that  the  zero  db  line  in  Fig.  19  lies 
below  the  flat  amplitude  ratio  portion  running  from  zero  to  the  breakpoint  at 
|l/Ts|,  this  increment  is  sufficient  to  make 
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thereby  reversing  the  original  assumed  inequality  (|Lp|  >  |LpN^/Np|)  responsible 
for  the  diverging  spiral  mode. 

Another  feature  of  a  pure  gain  yaw  rate  to  rudder  controller  is  its 
tendency  to  uncoordinate  turns.  This  can  be  seen  by  considering  the  side 
acceleration  equation  for  steady  state,  constant  altitude  conditions, 

ay  =  -g  sin  cp  +  U0r  =  Yyv  +  U0Y5*5r  (99) 

For  the  turn  to  be  coordinated  ay  should  be  zero,  so  that,  in  the  turn, 

Yvv  +  U0Yg*6r  must  also  be  zero.  When  ay  =  0,  yawing  velocity  and  bank:  angle 
are  then  connected  by  the  relationship 

r  =  —  sin  cp  (1  00) 

u0 

When  6r  =  0  and  Eq  100  holds,  Eq  99  implies  that  the  sideslip,  v,  is  also  zero. 
However,  when  8r  =  -Krr,  the  quantity  U0Yg*6r  will  have  a  value  other  than 
zero,  thereby  requiring  some  sideslip  if  ay  is  to  be  zero,  or  some  lateral 
acceleration  if  v  is  to  be  zero.  This  deficiency  of  a  pure  gain  r  -»«■  5r 
system  can  be  eliminated  either  by  introducing  an  rc  command,  or  by  using  pure 
lead— lag  equalization  ("washout")  in  the  r  &r  controller.  As  indicated  by 
the  Bode  and  root  locus  sketches  of  Fig.  20  the  introduction  of  lead— lag 
equalization  need  not  strongly  affect  the  augmented  Dutch  roll  damping, 
although  the  spiral  mode  is  no  longer  stabilized. 

2 .  cp  — ■-  5a  as  a  Single  Loop 

Single  sensor  control  loops  of  the  cp  8a  variety  have  been  thoroughly 
documented  elsewhere  (e.g.  Refs  2,  3,  and  13)-  Although  such  systems  can  exhibit 
a  wide  variety  of  behavior  depending  upon  the  relative  locations  of  the  poles  and 
zeros  of  and  the  equalization  used  in  the  controller,  the  most  usual  situa¬ 
tions  encountered  are  typified  by  the  family  of  root  loci  shown  in  Fig.  21  . 

Here  the  Dutch  roll  is  lightly  damped  and,  depending  upon  its  root ' s  location 
relative  to  the  £q,,  co^  numerator,  can  be  unstable  when  the  loop  is  closed.  With 
the  exception  of  these  "ou^,  a>d  effects",  the  single  axis  system  will  exhibit 


TO 


excellent  closed-loop  characteristics  when  simple  lead  equalization  (TEs  +  l) 
is  added  to  the  basic  airframe  characteristics.  This  equalization  can  be 
obtained  either  by  operating  on  the  cp  sensor  signal,  or  by  adding  a  roll  rate 
gyro.  In  the  latter  case  Tg  is  the  ratio  of  roll  rate  to  roll  position  gains. 

Several  means  (see  Ref  3)  can  be  employed  to  offset  deleterious  a>d 
effects  such  as  those  exhibited  when  the  $ga  numerator  quadratic  is  represented 
by  the  location  marked  with  the  "a"  subscript.  The  method  most  pertinent  to 
the  present  discussion  is  to  use  the  yaw  rate  loop  as  a  subsidiary  loop  to  the 
cp  — »-Sa  loop.  One  helpful  consequence  of  the  yaw  rate  closure  has  already  been 
noted  in  connection  with  Eq  96,  i.e.  the  increase  in  Dutch  roll  damping .  This 
will  result  in  the  Dutch  roll  poles  being  moved  further  into  the  left  half 


plane.  Another  beneficial  consequence  of  the  yaw  rate  closure  can  be  seen  by 


Figure  21  .  Root  Loci  for  Single  Axis  Bank  Angle  Control  System  with 
Ideal  Controller  and  Variable  Kqjg  >  Sa  =  K<p(TgS  +  1  ) 
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Using  for  Np  in  Eq  101  indicates  that  the  yaw  rate  closure  also  moves  the 

<^p  quadratic  further  into  the  left  half  plane.  In  the  roll  loop  closure 
the  locus  will  still  go  from  to  as  gain  is  increased,  hut  this  can  no 
longer  result  in  an  instability  if  the  roots  corresponding  to  ofy  and  0)d  are 
sufficiently  damped.  Since  this  is  accomplished  by  the  yaw  rate  loop,  it  is 
seen  to  be  an  admirable  inner  loop  for  the  cp  — < “-5a  system. 

C.  LOOP  CLOSURE  SEQUENCE 

With  the  benefit  of  the  discussions  above  a  unique  closure  sequence  can 
now  be  established  for  the  cpc  Ba,  r  5r  system.  Table  III  specializes 
the  general  factors  listed  in  Table  I  to  this  system.  The  net  results  of  the 
factors  summarized  in  Table  III  include  the  equivalent  block  diagram  of  Fig.  22, 
to  indicate  the  closure  sequence;  the  implicit  requirement  (which  goes 
alqng  with  the  absence  of  an  rc  in  Fig.  22)  for  washout  (i.e.  pure  lead— lag) 
equalization  in  Gr;  the  explicit  requirement  for  the  r  — ' — 5r  loop  to  appreciably 
augment  the  Dutch  roll  damping;  and  possible  lead  equalization  in  G^,. 

D.  NUMERICAL  EXAMPLE 


The  discussion  above  has  covered  the  general  aspects  of  the  bank  attitude 
control  system  with  a  washed-out  yaw  rate  inner  loop .  For  the  numerical  example 
of  the  present  article  a  second-order  lag  will  be  used  to  approximate  the  servo 
characteristics,  i.e. 
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2^5s 

-f  -  +  1 


(102) 


where  =  20  and  =  0.7  are  assumed  as  nominal  values.  The  same  servo 
characteristic  will  be  used  for  the  controller  transfer  function  in  both  roll 
and  yaw  axes . 

The  airframe  transfer  functions  are  slightly  modified  from  those  of  Ref  11 
to  provide  a  0§a  transfer  function  in  which  is  greater  than  oo^.  Also,  for 
this  example,  a  coupling  numerator  is  needed.  The  airframe  characteristics  to 
be  used  are. 
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Figure  22.  Equivalent  Block  Diagram  for  Bank  Angle  Control  System 
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Besides  the  vehicle  and  servo  characteristics,  values  are  needed  for  the 
yaw  loop  washout  time  constant,  Tp  and  the  roll  loop  lead-equalization 
time  constant,  T-g ,  Both  these  time  constants  would  ordinarily  he  found  as 
part  of  a  detailed  optimization  procedure  in  which  the  vehicle  dynamics  over 
the  entire  flight  regime,  and  the  effects  of  subsequent  loop  closures  (e.g. 
a  heading  loop),  are  considered.  Since  studies  of  this  kind  are  beyond  the 
present  scope,  reasonable  representative  values  are  selected  for  the  numerical 
example  .  With  these  values  the  controller  transfer  functions  become 


Gr  = 


Krs 


V~fo  +  1  ) 


Jcp 


17jL)2  g(.7)s  1 

|_20^  (20)  J 
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The  principal  objects  of  the  analysis  are  to  find  values  for  K^,  and  K® 
which  provide  ’’good”  closed-loop  characteristics  ("good"  implies  negligible 
static  errors  in  cp,  rapid  and  well  damped  responses  to  cpc,  adequate  margins  for 
vehicle  and  controller  pole  and  zero  variations,  relative  Insensitivity  to 
parasitic  nonlinearities,  etc.),  and  to  find  the  poles  and  zeros  of  the  closed- 
loop  transfer  function. 

The  first  step  in  the  analysis  procedure  is  the  simultaneous  closure  of 
the  yaw  rate  and  coupling  loops  to  obtain  $5a-  Bode  diagrams  for  Gr  /Ncpga 

and  GrRgr  are  shown  in  Fig.  2Ja  and  Fig.  2Jb  respectively.  The  closed-loop 
poles  resulting  from  closure  of  the  first  ( coupling  numerator)  will  be  the 
zeros  of  ®g'a,  while  the  closed-loop  poles  from  the  second  become  the  poles 
of  ®5a.  Only  Kr  is  available  for  adjustment  in  either  closure. 

Suitable  closure  criteria  can  be  expressed  in  several  different  ways,  e.g. 
as  phase  margins  and  gain  margins,  closed-loop  root  values,  integral  performance 
measures  such  as  ITAE  or  IE^  based  on  dominant  modes  (Ref  l6)  etc.  The  actual 
criteria  selected  are  somewhat  a  matter  of  taste  and  convenience,  but,  whatever 
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they  are,  the  resulting  system  should  possess  good  closed- loop  characteristics. 
For  the  present  case  the  closure  selected  will  be  based  on  design  to  a  specified 
damping  ratio  of  0.8  for  the  closed-loop  Dutch  roll  mode.  The  selection  of  a 
specified  damping  ratio  as  the  closure  criterion  provides  an  alternative  example 
to  the  use  of  a  phase  margin  criterion,  as  in  the  altitude  control  case.  The 
value  of  =  0.8  might,  at  first  glance,  appear  rather  high.  However,  the 
corresponding  value  of  obtained  from  the  simultaneous  closure  of  the  coupling 
loop,  will  be  somewhat  lower  than  0.8,  and  the  final  cp  ~*“"8a  closure  will  tend 
to  force  the  modified  Dutch  roll  roots  to  values  near  £  '  .  An  alternative,  of 
course,  is  to  specify  a  damping  ratio,  £',  for  the  3>'g  numerator.  As  it  will 

T  3. 

turn  out,  a  specification  of  =  0.8  is  equivalent  to  0.5  <  <  0.6. 

The  value  of  gain,  K^.,  consistent  with  a  of  0.8  can  be  found  in  one  of 
several  ways ,  e .g .  from  a  detailed  root  locus  plot .  In  the  present  example 
it  is  expedient  to  use  a  short  segment  of  a  so-called  g  Bode  plot  (see  Ref  9), 
for  £  =  -0.8,  as  shown  in  Fig.  25b.  The  value  of  |s|  where  4  G(g,  |s|)  is 
-l80  degrees  is  then  the  undamped  natural  frequency,  cd^,  of  the  Dutch  roll 
as  modified  by  the  yaw  rate  closure,  and  the  zero  db  line  location  is  found 
by  running  it  through  |G(|,  I s I ) | ^  at  this  value  of  | s | .  Kr  is  thus  found 
to  be  -8.88  in  linear  units.  (Kr  is  negative  because  the  sign  convention  for 
Sr  results  in  Nsr  being  negative.)  The  other  closed- loop  $5a  numerator  and 
denominator  roots  are  found  on  the  plots  of  Fig.  25  using  the  decomposition 
method.  The  resulting  $5a  transfer  function  is 


s'  -  JL 


8a~8a 


q 

a]. 


(-1 5,600)  (-§^  +  1) 

t,/ 

2(.58)s  „"l  T/  s  y2  2(.7)s  .1 

+  -b^TT  +  ll  +  TT8T2T  +  U 

(  -  .001 05  +  1  ^  .51  +  1  ^1 

.777  1 

,  s  n2  2(  .8)s  1 

_[l  .82^  (l  .82)  J 

fc  s.  ,2  2(.6?)s  1 

(_  18.5  (18.5)  J 

(105) 


The  last  quadratic  in  the  numerator  and  denominator  of  Eq  1  05  derive  from 
the  servo  characteristics .  These  are  nearly  identical,  and  can  be  presumed  to 
cancel  with  only  negligible  error.  This  type  of  near  cancellation  invariably 
occurs  when  coupling  numerators  are  present.  It  will  be  recalled  that  a  similar 
numerator  quadratic  was  not  present  in  the  altitude  controller  case  because  of 
the  absence  of  a  coupling  numerator  closure.  Cancelling  these  quadratics  gives 
the  result. 
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Diagrams  for  Bank  Angle  Control  System  Coupling  and  Yaw  Rate  Loops; 
GrNg  5  /%£  and  0^5 r,  Respectively 


®iL  = 


TP2“d2 


]h 


2 1  's 


<P 


+  — f-  +  1 

“<p 


] 


(^s  *  iKTJ,.  +  l)(TAs  *  .)|Q|)2  *  -|t-+  l| 


(106) 


-.00105 


+  D(~+1)( 


1  .777 


+  1 ) 


r.,2.  2(.8)s  .  . 

I  1  -|  Qq)  t\  Or\  \  +  1 


(1 .82) 
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The  ®8a  transfer  function  given  by  Eq  106  is  now  combined  with  (Eq  104), 

to  form  the  outer,  open-loop  transfer  function.  Bode  diagrams  and  a  root  locus 

sketch  for  closure  about  this  open-loop  transfer  function  are  shown  in  Figs.  24 

and  25.  To  make  the  closed- loop  transfer  function  cp/m  1 

cJ9c  —  5a 

r  — —  5r 

approximately  invariant  with  flight  condition  a  high  gain  closure  shall  be 
used  and  the  equalization,  (T£s  +  l),  will  actually  be  placed  in  the  feedback 

path,  i.e.  as  shown  in  Fig.  26.  The  closed-loop  transfer  function  will  then 
have  the  form 


JL 

% 


k- 

r 


(  v»«  \ 

T  *  Vi.  f 


Bq  (^TTJ  }  1  +  G^og 


(107) 


To  the  extent  that  the  high  gain  closure  makes  the  quotient  in  braces  approach 
unity,  tApJJ cpc  5a  "i11  approach  (tes  +  l)_1  .  T£  is  a  function  only  of 

r  —  5r 

controller  parameters,  so  its  value  can  be  controlled  readily. 

To  accomplish  the  high  gain  closure,  but  still  preserve  some  stability  (in 
the  closed-loop  servo  mode),  a  phase  margin  of  40  degrees  is  selected.  The 
gam,  K^,  will  then  be  -6.6  db  or  0.467  in  linear  units.  Using  unified  servo- 
analysis  techniques,  the  closed-loop  transfer  function  is 
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1  +  S°8a 


(l ) (TpgS  +  1)(Tes  + 

1) 

v  *  ^ 

+  1 

(T"s+ l)(T^s  +  l)(T^s+  1) 

(- 

c 

2  2?as 

i  +vr + 

a  ,2  -1 

[y  + 

(108) 


+  1  > <T§o  *  1  > 

fc  s  \2  2(  . 

‘-2.12  +~T2 

58)s  I 

TiiT  +  Ll 

.225  +  1  ^  .78  +  1  VlS.7  +  1  ^ 

(  s  Y 

_  2 .01  ' 

?  2(  -56)s  .  D  T/  s  \2 

+  723rr+1J[V.3) 

1 2(  -25)s ,  n 

TT^37  J 

Note  that  the  prime  notation  is  here  extended  to  two  loop  closures  by  the  double 
prime  quantities.  The  entire  numerator  is  approximately  cancelled  by  correspond 
ing  denominator  terms,  leaving 


1  + 


Vi 


'21  .7 


1) 


_  1 6 . ; 


2(-2l)s 

06.3) 


♦’] 


(109) 


This  is  nearly  unity  at  frequencies  below  12  rad/sec,  or  so,  and  hence 
^=]cpc  -*■-  8&  very  close  to  just  (TgS  +  l)  for  these  same  low  frequencies, 
r  — 5 

r 

The  width  of  this  frequency  region,  and  the  near  cancellation  indicated  by 
Eq  108,  cannot  be  materially  improved  without  serious  degradation  of  the  damping 
ratio  of  the  closed-loop  servo  mode. 

While  the  cp  response  to  <pc  signals  will  approximate  that  of  a  first  order 
system  having  a  time  constant  of  T£  =  l/o.8,  the  aircraft  response  to  external 
disturbances  will  not  be  of  such  simple  form.  For  example,  the  transfer  func¬ 
tion  relating  yawing  velocity  to  a  side  gust  will  not  contain  the  numerator  (to') 
quadratic  which  essentially  cancels  the  closed-loop  Dutch  roll  (a£)  mode,  as  ^ 
shown  in  Eq  108.  Instead,  the  dominant  yawing  velocity  mode  will  in  fact  be  the 
modified  Dutch  roll;  but  it  will  be  adequately  damped  (£”  =  .56).  Similar 
remarks  can  be  made  about  the  other  augmented-airframe  characteristics  which, 
in  general,  have  values  compatible  with  good  suppression  of  external  disturbances 
of  almost  any  type  likely  to  be  encountered. 
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Figure  24.  Bode  Diagram  for  Example  Bank  Angle  Control  System  Outer  Loop, 
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APPENDIX  A 


DEVELOPMENT  OF  EQUATION  24,  A,ys 

The  characteristic  function  of  the  generalized  multiloop  system,  is 

given  in  both  determinant  and  expanded  form  by  Eq  2k.  The  expansion  of  the 
determinant  is  considerably  expedited  by  writing  it  as  shown  in  Eq  A-1 ,  in  which 
each  of  the  element  columns  is  denoted  by  a  Roman  numeral  and  each  subcolumn  by 
a  capital  letter.  The  complete  expansion  of  Eq  A-1  produces  3^  =  8l  fourth- 
order  determinants  in  which  all  elements  are  single  (non-summed)  quantities. 
However,  as  will  shortly  be  seen,  60  of  the  determinants  are  zero,  one  is  equal 
to  Eq  25,  and  the  remaining  20  are  in  the  form  of  Eq  2 6  or  27. 


a  b  c 

a1  1  +  F11G11  +  F12G21 

a2l  +  f21g11  +  f22g21 

a31  +  F51G11  +  F32G21 

-a^i  +0+0 


II 

A  B  C 

a1  2  +  F1 1 G1  2  +  F12G22 
a~22  +  F2l  G1 2  +  F22G22 

a52  +  f51G12  +  f52g22 
•a42  +0  +  0 


III 

A  B  C 

al  3  +  F1  1  °1  5  +  F1  2°23 
agj  +  F21G,5  +  F22G2j 
ajj  +  F31G,j  +  FjgCfej 
-aj ^  +0  +  0 


IV 

A  B  C 

0  +  F-| ,  G,  ij.  +  F1  2G2^ 
0  +  Fg-]  G-|  +  F22g24 

0  +  F51Gl4  +  F5?G24 
1+0  +  0 


(A-1) 


Before  proceeding  with  the  systematic  expansion  of  Eq  A-1 ,  and  by  way  of 
review,  six  fundamental  properties  of  determinants  are  stated  (Ref.  12): 

1  .  If  all  elements  in  a  row  or  column  are  zero,  the  determinant 
is  zero. 

2.  If  all  elements  but  one  in  a  row  or  column  are  zero,  the 
determinant  is  the  product  of  that  element  and  its  cofactor. 

3.  If  two  columns  or  two  rows  of  a  determinant  are  identical, 
the  determinant  is  zero. 

4.  If  all  elements  in  any  column  are  multiplied  by  a  factor, 

•  the  determinant  is  multiplied  by  that  factor. 

5.  If  each  element  in  any  column  or  row  is  expressed  as  the  sum 

of  two  quantities,  the  determinant  can  be  expressed  as  the 
sum  of  two  determinants  of  the  same  order. 
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6.  If  two  columns  or  two  rows  of  a  determinant  are  interchanged, 
the  determinant  is  not  changed  in  absolute  value  but  is 
changed  in  sign. 

The  systematic  expansion  of  Eq  A-1  is  easily  accomplished  by  covering  two  of 
the  three  subcolumns  in  each  column  with  paper  strips  and  evaluating  the  visible 
fourth-order  determinant.  By  systematically  shifting  strips,  but  always  keeping 
two  subcolumns  in  each  column  covered,  which  is  equivalent  to  applying  property  5, 
Eq  A-1  is  completely  expanded  as  a  sum  of  determinants.  As  examples  of  the  evalua¬ 
tion  of  each  determinant : 


With  subcolumns  IA,  IIA,  IIIA,  and  IVA  exposed,  the  visible  determinant 
is  equal  to  A  (Eq  25). 


Shifting  one  strip  in  the  fourth  column  exposes  subcolumns  IA,  IIA,  IIIA, 
and  IVB,  and  the  visible  determinant,  using  property  4  and  the  definition 
of  NQj_  ,  is  G-itNq^  ;  i.e. 


+6i' 


a11 

a1 2 

a1  3 

f11g14 

a21 

a22 

a23 

f21  g1  4 

a31 

a32 

a33 

f31  g1  4 

'a4l 

_a42 

*a43 

0 

a11 

a1  2 

a1  3 

F11 

a21 

a22 

a23 

f21 

a  31 

a  32 

a33 

F31 

“a4i 

-a42 

_a43 

0 

g^HBi 


Shifting  one  strip  in  the  first  column  exposes  subcolumns  IB,  IIA,  IIIA, 
and  IVB,  and  the  visible  determinant,  using  property  4  and  then  property  5 , 
is  zero;  i.e.. 


F1 1 G1 1 

a1  2 

a1  3 

F1 1 G1 4 

F11 

a1  2 

a1  3 

F11 

F21  G1 1 
F31 G1 1 

a22 

a  32 

a23 

a33 

F21 G1  4 
f31  g1  4 

=  Gt  •,  G-,  4 

F2! 

f31 

a22 

a32 

a23 

a33 

F21 

f31 

0 

_a42 

"a43 

0 

0 

”a42 

-a43 

0 

In  general,  whenever  two  or  more  exposed  subcolumns  have  similar  alphabetical 
designations  other  than  "A",  the  determinant  is  identically  zero  by  properties 
5  and  4. 
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The  complete  set  of  possible  combinations  of  subcolumns  is  evaluated  in 

Table  A-I.  The  dots  in  a  particular  row  indicate  the  subcolumns  of  Columns  II, 

III,  and  IV  which  are  combined  with  each  of  the  three  subcolumns  of  I.  The 

evaluation  (nonzero  or  zero)  of  the  resultant  determinants  is  indicated  in  the 

appropriate  subcolumn  of  I.  The  nonzero  determinants  are  given  in  Table  A-II 

using  vehicle  transfer  function  numerator,  Nq<  ,  and  coupling  numerator, 
a-i  <1=  -  -  -  s.i 


notation  (Eq  2 6  and  27,  respectively), 


The  characteristic  function,  AgyS,  is  the  sum  of  all  determinants  shown  in 
Table  A-II,  and  is  given  below. 


of> 

w 

II 

> 

+  G11%,S1 

+ 

Gi 2Nq2& 

+ 

Gi  3r,q3gi 

+ 

G1  4wq45^ 

+  G2lNqi52 

+ 

g22%2&2 

+ 

G23Nq3g2 

+ 

G24wq4&2 

+ 

G1l( 

+ 

q1q2 

G22N51&2 

+ 

qlG3 

G23\4 

+ 

q4  \ 

Ga*VJ 

+ 

G12(G21\4 

+ 

<52q3 

+ 

) 

+ 

/  q3q1 

G13(G21N5fS2 

+ 

q3q2 

G22NBl62 

+ 

+ 

/  q4qi 

Gl4(G21Ngi52 

+ 

q4q2 

G22%52 

+ 

q4q3 

G25%62 

)  (A-2) 

Noting  from  property  6  of  determinants  that  NgigJ  =  and  recognizing  that 

controller  functions  other  than  those  chosen  can  be  factored  out,  it  is  clear 
that  various  simplifications  and  modifications  to  Eq  A-2  are  possible.  But  the 
simplest  of  all  is  the  shorthand  form  given  in  the  text  (Eq  24): 


4  2 

^sys  =  A  +  £  £  GjiKq 

1=1  j=1 


1BJ 


4  4 

£  £ 

i=1  k=1 
ijtk 


qiqk 

G1iG2kN5l52 


(A-5) 
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TABLE  A- I 


EVALUATION  OF  ALL  POSSIBLE  SUBCOLUMN  COMBINATIONS  IN  EQUATION  A-1 


I 

1 - m— 

IV 

A 

B 

c 

A 

B 

c 

A 

B 

c 

A 

B 

c 

S' 

S' 

S' 

• 

. 

. 

S' 

0 

s' 

■ 

. 

■ 

IS' 

S' 

0 

■ 

. 

. 

IS- 

0 

S' 

• 

. 

. 

0 

0 

0 

• 

• 

■ 

IS" 

0 

0 

■ 

. 

• 

S' 

S' 

0 

. 

. 

IS* 

0 

0 

■ 

. 

. 

0 

0 

0 

• 

• 

. 

S' 

0 

s' 

• 

• 

. 

0 

0 

0 

• 

. 

v*" 

0 

0 

. 

. 

0 

0 

0 

- 

• 

. 

0 

0 

0 

• 

. 

0 

0 

0 

• 

• 

0 

0 

• 

0 

0 

0 

iiBSIJ 

• 

0 

0 

0 

• 

• 

S' 

s' 

0 

- 

. 

S' 

0 

0 

• 

. 

0 

0 

0 

• 

• 

s' 

0 

0 

• 

. 

0 

0 

0 

■ 

• 

0 

0 

0 

B8MI 

■ 

. 

0 

0 

0 

PpS^j 

. 

. 

0 

0 

0 

■ 

- 

0 

0 

0 

• 

■ 

^  Denotes  that  determinant  of*  indicated  combination  of  subcolumns  is  nonzero 
0  Denotes  that  determinant  of  indicated  combination  of  subcolumns  is  zero 
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NONZERO 


EXPOSED  SUBCOLUMNS 
IA,  IIA,  IIIA,  IVA 
"  "  "  IVB 


"  »  "  IVC 


IIIB,  IVA 


"  "  "  IVC 


me,  iva 


"  "  »  IVB 


IA,  I IB,  IIIA,  IVA 


"  IVC 


"  IIIC,  IVA 


IA,  IIC,  IIIA,  IVA 


"  "  »  IVB 


"  IIIB,  IVA 


IB,  IIA,  IIIA,  IVA 


IVC 


IIIC,  IVA 


IB,  IIC,  IIIA,  IVA 


IC,  IIA,  IIIA,  IVA 


"  "  "  IVB 


IIIB,  IVA 


IC,  IIB,  IIIA,  IVA 


TABLE  A-II 


'S  IN  EXPANSION  OF  EQUATION  A-1 


VISIBLE  DETERMINANT 

A 

Gl  4IS5l 

G£*X2 

Gl  ^ 

Gi  3G2h\b2 

G^G23MB1&25 

G1  2%^ 

G1 2g24k51  B2 

G1 2g2JnB1  B2 

G22Nq2 

b2 

n  .  r  N%q2 

4g22nB1  B2 

gi3g22n6;52 

G1  1 

51 

q4 

G1  1  g24nb  b 

G1  °2 

q1  q3 

GnG25»BlB2 

q'.q2 

G1 1 g22^o  p~ 

G1  b2 

G2l  Nq, 

b2 

q4qi 

Gl4G2l\s2 

qBq1 

3°21  %  s2 

q2q1 

G1  2g21  5^ 
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APPENDIX  B 


TYPICAL  COUPLING  NUMERATORS 


Examples  of  longitudinal  and  lateral  coupling  numerators  are  given  in 
Tables  B-I  and  B-II,  respectively,  for  the  approximate  equations  of  motion 
shown  below: 

Longitudinal 


"s  -  Xu 

-Xw 

g 

0“ 

u 

Xg8 

-Zu 

( s  -  Zw) 

-U0s 

0 

w 

Z65 

-Mu 

-M^s  -  Mw 

2 

s  -  MqS 

0 

6  1 

= 

MgS 

0 

i/s 

-u0/s 

1 

h 

0 

Lateral 


- 1 

CO 

1 

K 

< 

-e/u0 

1  0" 

YgB 

-If3 

s(s  -  1^,) 

-L^  0 

9 

= 

I^B 

-Np 

-Nps 

s  -  N^.  0 

r 

NgB 

-UqS 

g  -  1XS2  sin  a0 

-UQ  -  lxs  COS  CLq  1 

ay 

0 

The  characteristic  determinants  of  both  the  longitudinal  and  lateral 
equations  of  motion  are  the  usual  quartics,  A^ong  and  with  unity  coeffi- 

1l 

cients  of  the  s  terms.  The  coupling  numerators  are  thus  compatible  with  the 
conventional  airframe  transfer  functions  in  Ref  2. 

In  the  tables  the  A,  B,  C  coefficients  are  those  for  terms  of  descending 
order  in  s  (e.g.,  A  is  the  coefficient  of  either  s2  or  s ;  B,  the  coefficient  of 
s  or  s°;  C,  the  coefficient  of  s°,  or  nonexistent).  The  complete  literal 
expressions  for  each  coefficient  are  shown,  as  are  the  forms  of  the  coupling 
numerators,  and  approximate  literal  values  of  the  gains,  time  constants,  fre¬ 
quencies,  etc.,  appropriate  to  these  forms.  In  most  cases  the  approximations 
involved  are  obvious  by  reference  to  the  complete  expressions  and  stem  largely 
from  the  usually  negligible  nature  of  certain  control  effectiveness  derivatives 
(e.g.,  Xg  longitudinally  and  laterally)  and  products  of  cross-control 
effectiveness  derivatives  (e.g.,  NgaLgr  «  LgaNgr)  . 
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One  property  of  coupling  numerators,  i.e.,  =  -]$£*  is  used  extensively 

in  the  tables  to  simplify  the  presentation.  An  additional  property  applicable 
for  the  determination  of  coupling  numerators  containing  an  auxiliary  output 
motion  quantity  is  also  used.  With  q.  (i  =  1,2, 3)  the  independent  degrees  of 
freedom  and  g^  an  auxiliary  variable. 


then 


%  aipi  qi  +  a4 2q2  +  a43q3 


„'Vi  ,  ..*1% 


'  s*i'»S -  v£<2  * 


(B-5) 


WaSae  °ne  °f  thS  couPllnS  numerators  on  the  right  side  will  be 


N8k5i  1S  zero  if  i  =  j  or  k  =  1 


zero  because 


For  example,  the  independent  degrees  of  freedom  for  the  longitudinal 
equations  of  motion  are  u,  w,  and  0,  and  two  auxiliary  motion  quantities 


are 


Using  Eq  B-3  and  B-4, 


h  = 


Uofl  1 
~9~  sw 


a 

6e6T 


—  N0  w 
UD  W6e8T 


Similarly,  using  Eq  B-3  and  B-5, 


N 


9  h 

5>e®T 


9  w 
5e8T 


Thus  from  Eq  B-6  and  B-7, 


sN 


0  h 

Be6j 


-UqN 


0  a 
Se&T 


(B-4) 


(b-5) 


(b-6) 


(b-7) 


(B-8) 
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The  property  of  Eq  B-3  that  coupling  numerators  can  be  summed,  to  yield  new 
numerators  applies  even  when  an  auxiliary  output  motion  quantity  is  expressed 
in  terms  of  other  auxiliary  variables.  If  Eq  B-5  is  written 


then 


or 


h 


a 


£ 

sesT 


IT  ^BeST 


0  a 
5eBT 


slJ, 


.a  h 
SeSrp 


TT  CL 

°  6e5T 


Comparison  of  Eq  B-8  and  B-11  shows  that 
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TABLE  B-II 

LATERAL  COUPLING  NUMERATORS 
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APPENDIX  C 


AN  ALTERNATE  FORMULATION 
OF  NUMERATOR  AND  DENOMINATOR  EXPANSIONS* 

Although  the  mathematical  equations  for  the  effective  (closed-loop)  system 
transfer  function  numerators  and  denominator  may  appear  quite  complex,  they 
are  simple  straightforward  expansions.  These  expansions  can  be  defined  by  a 
few  rules  similar  to  the  ones  used  in  signal  flow  diagrams.  The  sensing, 
actuation,  and  equalization  components  of  the  system  are  treated  as  feedbacks, 
all  acting  to  modify  the  vehicle  open-loop  transfer  functions.  The  resultant 
closed-loop  system  is  analogous  to  an  effective  vehicle  which  may  have  command 
control  deflection  or  disturbance  inputs. 

The  rules  are: 

1.  The  effective  numerator  is  equal  to: 

a.  The  open-loop  numerator 

b.  Plus  the  sum  of  all  the  feedback  transfer 
functions,  each  one  multiplied  by  the 
appropriate  coupling  numerator 

2.  The  effective  denominator  is  equal  to: 

a.  The  open-loop  denominator 

b.  Plus  the  sum  of  all  the  feedback  transfer 
functions,  each  one  multiplied  by  the 
appropriate  numerator 

c.  Plus  the  sum  of  all  the  feedback  transfer 
functions  taken  two  at  a  time,  each  pair 
multiplied  by  the  appropriate  coupling 
numerator 

Rule  1  simply  states  that  each  feedback  modifies  a  numerator  by  adding  to 
that  numerator  the  product  of  the  feedback  transfer  function  and  the  appro¬ 
priate  coupling  numerator.  The  appropriate  coupling  numerator  is  the  one  for 
the  output/input  pair  of  the  original  numerator  plus  the  feedback  output/input 

pair.  For  example,  the  feedback  qp— < *-6.  modifies  the  q-/&i  numerator  by  adding 

q .  q?  J 

to  it  the  term  G.  an<^  modifies  the  q. /r).  numerator  by  adding  to  it  the 

q .  qD  J  '  1  J 

term 

^Developed  by  Robert  L.  Stapleford. 
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The  effective  denominator  has  been  shown  in  Appendix  A  to  be  given  by 

4  2  4  4  q.q, 

^sys  =  A  +  Z  Z  Gji%iR  +  E  E  GliGa^Sp  (C-D 

i=1  j=1  i=i  k=1  1  2 

The  equivalence  of  Eq  C-1  and  Rule  2  should  be  obvious.  The  appropriate 

numerator.  Rule  2-b,  is  the  one  with  the  same  output/input  pair  as  the  feed¬ 
back  and  the  appropriate  coupling  numerator.  Rule  2-c,  is  the  one  with  the 
same  two  output/input  pairs  as  the  feedbacks. 

For  example,  consider  the  system  shown  in  Fig.  C-1,  which  is  the  same  as 
in  Fig.  4,  redrawn  to  form  the  command  control  deflection  input,  BP  (the 

‘-Q. 

product  of  the  command  input,  q-^  ,  and  the  command  controller  feedback 
transfer  function,  Gqi )  • 


Figure  C-1  .  Block  Diagram  of  Multiloop  System 
l|c  ’  'llj.— *-52'  "~B1 
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can  be  written 


The  closed- loop  system  transfer  xuin-c'.c 


q1/q1c 


iL 

qic 


<ii 


J21 


82 

FiT— 


(C-2) 


The  effective  numerator  for  the  system  with  four  loops  closed  is  determined  by 
applying  Rule  1 . 

r  +  r 

G1*t-N8281 


"<ll 


82 


"l! 


&2 


(C-3) 


where  the  two  bracketed  terms  are  zero  because  is  zero  if  i  -  j  or  k  =  1 

Using  Rule  2, 


A'"'  -  Agys  =  A  +  Gi^  +Gi4%6i  +  +  C21% 

Gl2fJl4N8161J  +  G1  2c2I+N61  82 

r  ^ 

'  G- 2^21^8^82  +  1^1  4G24%i  02J 

qk<ii  r  vil 

+  Gi4g21n8i62  +  |^2^G211,8282j 


where  the  three  terms  in  brackets  are  zero  because  of  the  coupling  numerator 
property  stated  previously. 

While  the  above  discussion  and  example  and  the  derivations  in  the  body  of 
the  report  have  considered  only  the  cases  of  the  command  response  with  feed¬ 
backs  to  two  independent  contiols  or  Uie  disturbance  response  with  feedbacks 
to  one  control,  a  complete  generalization  (for  a  three  degrees  of  freedom 
system)  is  quite  straightforward  with  the  introduction  of  the  coupling-coupling 
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numerator.  The  coupling-coupling  numerator  has  three  columns  of  the  open-loop 
characteristic  determinant  replaced  by  the  appropriate  control  or  disturbance 
coefficients. 


frile  1  is  generalized  by  adding  "c.  Plus  the  sum  of  all  the  feedback 
transfer  functions  taken  two  at  a  time,  each  pair  multiplied  by  the  appropriate 
coupling— coupling  numerator* M 

Rule  £  is  generalized  by  adding  "d.  Plus  the  sum  of  all  feedback  transfer 
functions  taken  three  at  a  time,  each  combination  multiplied  by  the  appropriate 
coupling-coupling  numerator." 


Por  example,  with  feedbacks  q.^6,,  qg  —  6g,  and  qj— the  disturbance 
response,  q-,/n,  is  given  by: 


n 


(C-5) 


A***  = 


*  *  [“nsrfajej] 

A  +  S11  %  6i  ♦‘*^♦<•3*^ 

+  G! !  g22nb]  6g  +  G11G33^63  +  g22G33"6^6j 


(C-6) 


+  G 


-  „  ql  *2*3 
nG2263JN51S£~3 


(C-7) 


where  tne  bracketed  terms  of  Eq  o'  -..re  fete  vat  a  coupling-coupling 

numerator  is  zero  if  any  two  of  the  outputs  or  any  two  of  the  inputs  are 
identical. 
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